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We consider the relationship between the tight-binding Hamiltonian of the two- 
dimensional honeycomb lattice of carbon atoms with nearest neighbor hopping only 
and the 2 + 1 dimensional Hamiltonian of quantum electrodynamics which follows in 
the continuum limit. We pay particular attention to the symmetries of the free Dirac 
fermions including spatial inversion, time reversal, charge conjugation and chirality. We 
illustrate the power of such a mapping by considering the effect of the possible symme- 
try breaking which corresponds to the creation of a finite Dirac mass, on various optical 
properties. In particular, we consider the diagonal AC conductivity with emphasis on 
how the finite Dirac mass might manifest itself in experiment. The optical sum rules for 
the diagonal and Hall conductivities are discussed. 

Keywords: Graphene; Dirac theory; chirality; symmetry 



Preprint: NSF-KITP-07-126 

i 



February 1, _ 2008 11:55 WSPC/INSTRUCTION 

review-graphene'jversion 



FILE 



2 V.P. GUSYNIN, S.G. SHARAPOV, J. P. CARBOTTE 
1. Introduction 

A single sheet of carbon atoms tightly packed into a two-dimensional (2D) honey- 
comb lattice is called graphene. The band structure of graphene consists approx- 
imately of a valence (full) and conduction (empty) band both conical in shape 
with vertex meeting at a point called a Dirac point. There are two inequivalent 
pairs of such cones. This model, which is based on the continuum limit of tight- 
binding bands can adequately describe low energy phenomena. Deviations from 
simple cones can become important however near the ends of the bands where 
additional details of the tight-binding bands need to be considered. An example 
of the need to go beyond the continuum approximation, which will be considered 
in this review, is a discussion of optical sum rules which involve the integration 
of the optical conductivity or closely related quantity over all energies within the 
band. An important aspect of the graphene problem is that it can be mapped, in 
the continuum approximation, into the Hamiltonian of 2 + 1 dimensional quantum 
electrodynamics (QED2+1) with Dirac fermions. It is remarkable that historically 
this field theoretical view of graphene started 20 years before^ its actual discovery^ 
and has motiv ated theoretical work on a condcnscd-mattcr realization of the parity 
anomalypEHl the renormalization group approach in this systemf^El anc i npn per- 
turbative dynamics of the generation of an excitonic-like gap in graphene 1SI9J T ne 
interest in Dirac fermions in condensed matt er theory was also related to studies 

IQI 

m 

The experimental proof of the existence of Dirac fermions in graphene^ESl 
that came fr om the observation of the theoretically expected unconventional 
with the quantized filling factor 



of the integer quantum Hall eff ect (IQHE )ff^ Nowadays graphene inspires work on 
the fractionalization of f ermiona ^ 12 | 13 | w ] 1 i c ] 1 s \j\\ waits for the discovery. 



<rxy = -rv, i/ = ±2(2n+l), n = 0,l,... (1) 

has promoted Dirac fermions from a beautiful theoretical toy to a real object that 
one day may perform in a " graphenium inside" processor-^ or even s oone r become 
a resistance standard operational above liquid-nitrogen temperature. 1 ^ This for- 
tunate situation has created much excitement and it is now possible to consider 
doing bench top experiment to explore the properties QED 2 +i. This mapping into 
QED also means that insight obtained in the study of relativistic fermions can be 
brought to bare on the study of a solid state system. Examples of relativistic effects 
that could be studied include the Klein paradox w hich deals with the perfect trans- 
mission of electrons through high potential and Zitterbewegung or 
trembling of the center of a free wave packet!^ The experimental realization of 
graphene calls for an overview of the theoretical link between the condensed matter 
lattice description of graphene and its continuum QED2+1 formulation. This should 
help one judge better which theoretical concepts developed during previous theo- 
retical studies are important for a deeper understanding of the analogies between 
the behavior of electrons in graphene and in QED2+1. 
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In this short review we consider in detail this correspondence between a tight- 
binding lattice model and QED2+1. We concentrate on symmetries of free Dirac 
fermion problem and the consequences of possible lacking of such symmetries. In 
particular, we consider spatial inversion, time reversal, charge conjugation and a 
continuum U (4) symmetry. We explain the difference between 2D Dirac fermions in 
graphene and the 3D Dirac fermions that are studied in the more familiar QED3 + i. 
In particular, there is a difference in the physical meaning of such an important 
quantum number as chirality in QED2+1 and QED3 + i. Because of this and due to 
the fact that the notion of chirality is now widely used in the literature on graphene 
(see e.g. Ref . l25f . we found that it is useful to explain the similarities and differences 
of this notion in field theory and graphene. 

Another important issue is the possible effect of correlatio ns on various proper- 
ties of graphene. For example, the new IQHE statePETEHl 

with the filling factor 

v = 0,±1 seen in the strong magnetic field B greater than or equal to 20 T are 
likely to originate from many-body interactions. In partic ular, the zero filling fac- 
tor [y — 0) state is related to a spin polarized stat cF^ -^, while v = ±1 states are 
probably related to the lifting of sublattice degeneracy. There also exists evidence^ 
that Coulomb interactions may be partially responsible for an observed increase in 
the energy difference between second and first Landau levels in fields ~ 18 T. For 
the conventional case a theorem by KohrPSl based on nearly parabolic bands shows 
that there should be no such shifts due to correlations, but the theorem may not 
apply here. An interesting feature of the QED description of graphene is that there 
is a phenomenon of spontaneous symmetry breaking when the interactions generate 
a Dirac mass or gap. Historically this was exactly the interest in the condensed 
matter realization of the an omalous properties of massive QED2+1, which brought 
the attention of theoretician a^^'^ to graphene before its actual discovery. Now it 
is an open question (see e.g. Ref. [31] for a brief overview and more references are 
provided below in Sec. 13 . 6[) whether the spontaneous mass generation and symme- 
try breaking indeed occur in graphene or the observed v = ±1 IQHE states have a 
different origin. Because the experiment does not yet provide a definitive answer, in 
this review we discuss the theoretical possibility of breaking the U (4) symmetry of 
the Lagrangian of graphene by eight Dirac masses. We illustrate how the creation of 
a finite Dirac mass affects the discrete and continuum symmetries and trace in one 
case its consequence in the diagonal AC conductivity. It is left to future experiments 
to verify this theoretical possibility. 

Although graphene is only 3 years old, there already exi sts several popular ar- 
ticles at the introductor y^ ! ^ and more 

advancecP13513U 

level. Also there is an 

excellent review article by A.K. Geim and K. Novoselov^ and a more theoretical 
review by M.I. Katsnelson and K.S. NovoselovPS Because graphene is a basic ele- 
ment of all graphitic forms, all existing literature from late 40's onwards considers 
graphene as a starting point of the analysis. Particularly useful are the results for 
carbon nanotubes (see e.g. Refs. l36 | [37 ]) . 

Despite this variety of works on graphene we found that a more formal view 
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Fig. 1. (Colour online) (a) Graphene hexagonal lattice constructed as a superposition of two 
triangular lattices A and B, with basis vectors ai ^ for lattice A and vectors Si with i = 1,2,3 
connecting A to B. (b) The green hexagon is a Brillouin zone (BZ) and pink diamond is the 
extended BZ for the honeycomb lattice. The reciprocal lattice vectors are bi and b2. 

which can form the bridge between condensed matter and quantum field theory is 
still missing, and we decided to present it here. The review is organized as follows. 
In Sec. [21 we introduce the tight-binding model with an external vector potential and 
obtain expressions for electric current and diamagnetic term. In Sec. [3] we derive 
the continuum QED2+1 Lagrangian from the tight-binding model of graphene. This 
allows one to trace the link between the underlying lattice structure of graphene 
and clarify the physical meaning of the spinor components in the Dirac theory. The 
physical meaning of chirality in graphene is discussed in Sec. 13.41 In Sec. 13.51 the 
discrete symmetry operations of the effective Dirac theory are defined in accordance 
with the discrete symmetries of the lattice model considered in Sec. 12.31 Basing 
on the discrete symmetries we discuss in Sec. 13.5.41 the difference between massless 
neutrinos in QED 3+ i and quasiparticles in graphene. Eight possible Dirac masses are 
introduced in Sec. l3.6l and their transformation properties under discrete symmetries 
are classified. The AC conductivity in zero and finite magnetic field and optical sum 
rules are considered in Sec. [4] In Sec. [5] we present our conclusions and discuss open 
questions. 

2. Tight-binding description of graphene 
2.1. Tight-binding model 

The honeycomb lattice can be described in terms of two triangular sublattices, A 
and B (see Fig. [T]a). A unit cell contains two atoms, one of type A and one of type 
B. The vectors 

& i = a (^—)^ & 2 = a(-, — — ), (2) 

shown there are primitive translations, where the lattice constant a — \a.± \ — |&2 1 = 
\^3acc and ace is the distance between two nearest carbon atoms. The correspond- 
ing reciprocal lattice whose vectors are bi = ^£(1, 1/V3) and b2 = (1, —l/s/3) is 
shown in Fig.[T]b together with the reduced (symmetrical and extended) Brillouin 
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zone. The reciprocal vectors satisfy the relation • bj = 2^8^. 

Any A atom at the position n = aini + a2«2, where ni,n2 are integers, is 
connected to its nearest neighbors on B sites by the three vectors Si: 

Si = (a x - a 2 )/3, <5 2 = ai/3 + 2a 2 /3, <5 3 = -5i - 5 2 = -2a x /3 - a 2 /3. (3) 

Besides translations, the symmetry group of honeycomb lattice includes rotations (R 
and R^ 1 ) on ±27r/3 and mirror reflections (Yy, Y%,Y-i) about planes passing through 
the center and three hexagon vertices. Together the rotations and reflections form 
a nonabelian group C% v with 6 elements. Their explicit matrix representation may 
be defined by matrices giving the transformations of the basis vectors ai and a 2 , 

*-(- 1 1 0-*-Gi)- l4 -(-°.~o 1 )- 

In addition there is a reflection Z in the plane of the sheet (changes a^ — ► —a,;, 
followed by the translation by one of the vectors 38i) and, accordingly, combina- 
tions of Z with any of the above operations. All these operations do not interchange 
A- and B-type atoms. The one-electron eigenfunctions can be classified according 
to the subgroup (E, Z) whether the states are even (a states) or odd (tt states) 
under reflection Z. There are also operations which interchange A and B atoms, 
for example, reflections Xy,X2,X$ in mirror planes perpendicular to the corre- 
sponding Yi planes. They are symmetry operations when accompanied by some 
fractional translations!^ The rotation and rotation-reflection operations (leaving 
aside Z reflection) form the point-group of graphene which contains 12 elements. 
The effect of group operations G on a function of the coordinates is defined as 
ip'(r) = -0(G _1 r) = T(G)tp(r), and it allows one to determine the irreducible rep- 
resentations for graphene. In particular, one should note that the subgroup ((4]) has 
a two-dimensional spinor representation which we use below. Full consideration of 
the symmetry group of graphene was made in Refs. [33 l39l 

The carbon atoms in graphene plane are connected by strong covalent er-bonds 
due to the sp 2 hybridization of the atomic 2s,2p x ,2p y orbitals. The 2p z (n) or- 
bitals are perpendicular to the plane and have a weak overlap. Therefore we start 
with the simplest tight-binding description for n orbitals of carbon in terms of the 
Hamiltonian 



n,5i ,<7 



al <(7 exp ( r-<5jA ) b n+S:Cr + c.c. 



(5) 



where t is the nearest neighbor hopping parameter, a ni<T and 6 n +6,cr are the Fermi 
operators of electrons with spin a =t,| on A and B sublattices, respectively. 
Since we are interested in the current response, the vector potential A is intro- 
duced in the Hamiltonian (J5]) by means of the Peierls substitution <v n a b m . a — > 
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a' a a exp (— || f m Adr) 6 m>0 . that introduces the phase factor exp(||5jA) in the hop- 
ping term (see Ref . l40l for a review). We keep Planck constant h and the velocity of 
light c, but set kg = 1. The charge of electron is — e < 0. Note that we consider a 2D 
model. The experiment shows that graphene, in order to exist without a substrate, 
spreads itself out of the 2D plane!^ 

Expanding the Hamiltonian (|5]) to the second order in the vector potential, one 

has 



H = H Q -Y / 



~A(n)j(n) - ^A a (n)T al3 (n)Ap(n) 



a, /3 = 1,2. (6) 



The total current density operator is obtained by differentiating Eq. © with respect 
to A a (n) , 

dH 



j a (n) = -- 



d(A a /c) 

and consists of the usual paramagnetic part 



= 3 a (n)-T a 0(n)Ap(n)/c, 



(7) 



(8) 



and diamagnetic part, 

f ^ &2H 
T a /3{n) 



te 2 



d{A a /c)d{Ap/c) fi? 



^2(Si) a (Si)(3 al^bn+s^ + K+S,* ™," ■ (9) 



We stress that to obtain the correct form of the current ([8]) and the diamagnetic 
part © , the Peierls substitution has to be made in the initial Hamiltonian ([5]) or its 
momentum representation (see Eq. (|10| below) rather than after the diagonalization 
of these Hamiltonians is made (see Ref. and references therein). We will return 
to the paramagnetic and diamagnetic terms of the interacting Hamiltonian (|6|) later 
when we discuss coupling to the vector potential in Sec. 13.21 and the sum rules in 
Sec. 14.41 For now we consider the noninteracting Hamiltonian Hq. 



2.2. Spinor representation of the noninteracting Hamiltonian 

In the momentum representation the Hamiltonian Hq reads 

with 0(k) = -tJ2 Sl eik<5, = -e(k)e ll ^ k ). In the basis © we get 
0(k) = -te ik( - ai - a2)/3 [1 + e ika2 + e~ lkai ] 



exp i 



. kyQi 

V3 



exp 



_ , 2 cos — — 

2V3J 2 



(10) 



(11) 
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e/t 



-2 




♦ x v .. 




ak„ 



ak x 



Fig. 2. (Colour online) The energy band structure of graphcne. Valence and conduction bands 
meet at six K points. 



and 



e(k) = t\ 1 + 4 cos 2 h 4 cos cos 



(12) 



2 2 2 

Accordingly, because the graphene structure contains two atoms per unit cell (two 
sublattices), the s pect rum of quasiparticles excitations has two branches (bands) 
with the dispersiorJ^I E± = ±e(k) shown in Fig. [2] In Eq. (fT0)l we introduced the 
spinors 

'a<r(k)" 



& CT (k) 



(13) 



with the operator T CT (k) being the Fourier transform of the spinor T CT (n) = 



K.o 



T (7 (n) = VS [ ^e ikn T CT (k). 



(14) 



Here 5 = v / 3a 2 /2 is the area of a unit cell and the integration in Eqs. (fTU|) and 
(fl"4"l) goes over the extended rhombic Brillouin zone (BZ). We also add to Ho the 
Zeeman term and the chemical potential 

d 2 k 



Tt(k)T CT (k) 



(15) 
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with fie = fj, — where fig = eh/ (2m e c) is the Bohr magneton, and g is 

the Lande factor. We count the Zeeman energy from the chemical potential fi, so 
that our subsequent consideration is based on the grand canonical ensemble. 

The corresponding imaginary time Green's function (GF) is defined as a thermal 
average 

G a ( n -T 2 ,n 1 - n 2 ) = -(TrT^Ti.niyrtfo.na)) (16) 

and its Fourier transform is 

G a (n -T2>ni -n 2 ) = (17) 
d 2 k 



G CT (ia;„,k)exp[-zw„(ri - t 2 ) + ik(ni - n 2 )] 

BZ ( 27T ) 



with 



G,(z^,k) = ( ^ + ^ )2 _ £2(k ) u» = *(2n + l)T, (18) 

where the matrices t± = (tj ± iT 2 )/2 made from Pauli matrices, operate in the 
sublattice space. The GF (fT8"l) describes the electron- and hole-like excitations with 
the energies E±(k) — ±e(k) — u a , respectively. 

The dispersion e(k) near the six K points ±27r/o(l/3, l/y/S), ±27r/a(2/3, 0), 
±27r/a(l/3, -1/V3) at the corners of the hexagonal BZ (see Fig. [2]) is linearp^ 
£-±(p) = ±ft«F |p|— Mcr (see Fig.[3](a) in Sec. l3.6p . where the wave vector p — (pi,p 2 ) 
is now measured f rom th e K points and the Fermi velocity is vf = y/3ta/{2h). Its 
experimental value^IISl i s Vp ~ 10 6 m/s (see also Ref . Hll for more recent values of 
Vf and t). Since only two K points are inequivalent, in what follows, we select them 
to be K± = ±27r/a(2/3, 0) (see Fig. [T]b), so that they are inside the extended BZ. 
The degeneracy of the two K-t points is protected by the point-group symmetry of 
the hexagonal lattice. 



2.3. Discrete symmetries of the tight-binding Hamiltonian 

The K± points in single layer graphene are stable against perturbations which 
preser ve t he discrete spacetime inversion symmetry and which do not mix the K-t 
points3^1 We consider these symmetries for the tight-binding model and define 
below the corresponding discrete symmetry operations for the effective QED 2+ i 
theory of graphene. A topological stability for the appearance of massless Dirac 
fcrmions was studied in Ref. and in Sec. 12.3.31 we discuss it. 



2.3.1. The spatial inversion V 

Choosing the center of symmetry to be the center of the hexagonal in Fig. Q] a, we 
observe that the spatial inversion V: (x,y) — * (— x, —y) would be the symmetry of 
the system if one also exchanges A and B atoms, i.e. 

a n ,<7 -* Van^V^ 1 = 6_ n>CT & n>(T -> Vb^V^ 1 = a_„ !(T . (19) 
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As we will discuss in Sec. 13.51 the definition of V is different from that commonly 

11 Ql 

used in QED2+i. LL ^ J In momentum space V acts as follows 

a CT (k) - Va^V- 1 = 6 CT (-k), b a (k) - Pb.Qn)?- 1 = a CT (-k), (20) 

i.e. it reverses the sign of the momentum k — > k and exchanges K points. For the 
spinors T cr (k), the action of V is defined as 

T CT (k) — > T'T^k)^- 1 = nT CT (-k). (21) 

One can easily check that the noninteracting Hamiltonian Hq is invariant under V 

H — ► T'i/o^ 1 = #o (22) 

because the Hamiltonian density 7io satisfies the condition 

TiWo(k)Ti=Wo(-k). (23) 

However, if one assumes that the densities of particles on the A and B sublattices 
are different, i.e. that 

ff i = E/ 7^K4(kK(k)+m 6 &t(k)Mk)] (24) 

„ JBZ l Z7r J 



E 



/ 7T7^Tt(k)[m + r + m_r 3 ] T ff (k) 
Jbz 



with m + = (m a + m&) /2 [tq is 2 x 2 unit matrix] and m_ = (m a — m&) /2 such term 
would break the inversion symmetry by swapping now inequivalent sublattices: 

Hi — > Pi/^- 1 = V / K4(k)«a(k) + m a b\ (k)6 ff (k)] (25) 

or in spinor notation 

if! — = £ (k)(m + r - m_r 3 )T CT (k). (26) 

One can see that the term with m + that corresponds to the same carrier density on 
A and B sublattices can be absorbed in the chemical potential fi, and the carrier 
imbalance term with m_ is parity breaking. In practice, it has been suggested that 
one can introduce a Dirac mass term by placing the graphene sheet on a substrate 
made of hexagonal boron nitride. In this circumstance the two carbon sublattices 
become inequivalent because of its interaction with the substrate. A recent band 
structure calculation^ for such a configuration has given a gap of about 53 meV. 



2.3.2. Time reversal T 

The time-reversal operation, t — > — t changes the signs of the momentum and spin 

AO 

leaving the sign of coordinates unchanged!-* 21 

/a+(k)\ /a+(k)\ x ( a _(-k) \ /a+(-k)\ 

UwJ^Uwr =U + (-k)J= jCT2 U(-k)J' (27) 
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The operator b± (k) obeys the same rule, and a± (k) , b± (k) transform as follows 

(4(k) a!(k)) ->r(aj.(k) a^k))^ 1 = ( a J.(k) ol(k)) H<r 2 ). (28) 

To invert the direction of time, the operator T has to be antiunitary. The action of 
T on the sublattice spinors is given by 

T(k) — ► TT(k)T -1 = i<r 2 T(-k), (29) 

where oi acts on the spin indices of the spinor Y CT (k). 

One can check that the noninteracting Hamiltonian Hq is invariant under T 

Ho — > THqT^ 1 = H (30) 

because 7io satisfies the condition 

roW5(k)ro =W (-k). (31) 

Since a fixed external magnetic field breaks time reversal symmetry, both the Hamil- 
tonian ([SJ which includes this field and the Zeeman term (|15[) break it. On the other 
hand, the Hamiltonian Hi given by (|24[) is invariant under time reversal, but if we 
consider spin dependent masses m aj b !CT , it will also break this symmetry. 
The presence of combined VT symmetry enforces a condition on Ho : 

T 1 n* (k) n = n (k) (32) 

or in matrix components: Hq 1 = Hq 2 and H}, 2 = H 21 * ■ One can say that the first 
equality forbids m_T3 terms in the Hamiltonian protecting K± points when they 
are not mixedJ^l The energy spectrum is E = Hq 1 ± |^o 2 l- For H^ 1 — const, a 
constant term can be absorbed in the chemical potential, so that VT symmetry 
results in the symmetry of the spectrum similarly to the relation (|33p considered 
below. 

2.3.3. Particle-hole symmetry of the spectrum and stability of Fermi point 

The Hamiltonian Ho(k) also satisfies two other relationJ3^HESISni ( see a i so the 
end of Sec. I53JI 

T 3 H (k)T 3 = -Wo(k) (33) 

and 

r 2 7i*(k)r 2 =-W (k). (34) 

It is easy to check that they guarantee that if there is a state with energy 
E, then the states Ts\ip) and T2\4>}* correspond to the energy —E which implies 
that the energy bands are symmetric about E = 0. We note that while the parity 
breaking term m-T^T obviously violates the condition ([33]), the more general 
condition (|34[) is still satisfied. Indeed, even for to_ ^= the spectrum is symmetric 

about E = 0: E(k) = ±Jm 2 _ + \4>(k)\ 2 . In order for the Fermi point E(k = 0) = 
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11 



to exist, we must require the relation (|33p which is a necessary but not sufficient 
condition for its existence. The existence and stability of the Fermi point is dictated 
by topology in momentum space^U that is by nonzero topological invariant (winding 
number) which is expressed analytically as 

N= i ^M[a s H- l ()L)d±H{\i)). (35) 



Here the integral is taken over an arbitrary contour C around one of the K.± points 
and tr is the trace over sublattice indices. For the Hamiltonian Ho given by (|10ll 
the topological invariant can be rewritten in the form 

1 ^ 0(k)d0*(k) - <f (k)#(k) 



N 



(36) 



ic \m\ 2 

For K± points this gives 2V(K±) = ±1. On the other hand, when the spectrum 
becomes gapped, the topological invariant N = and the Fermi surfaces disappear. 
This can happen because the total topological charge at two K± points N = 0. The 
trivial total topological charge of the Fermi surfaces allows for their annihilation, 
which occurs when the energy spectrum becomes fully gapped. 

Now we come to the question of topological stability of the Dirac points with 
respect to the presence of the other hopping terms.SSl 

The next neighbor (intrasublattice) hopping considered, for example, in Rcf. 52, 
introduces only the diagonal term in the Hamiltonian 7io in Eq. (jlOj) which elim- 
inates the particle-hole symmetry of original energy bands and for which small 
deviations from K± can be absorbed in the chemical potential!^ Thus, we address 
here the role of the diagonal transfer (third neighbor hopping) t' which is described 

by 



H t > = -t'^a^bn^S! +c.c. 

Accordingly, in the presence of Hf, Eq. (fTTjl for (j>(k) takes the form 
<£(k) = -i^y k<5 * ~t'e- 2lkSl 



(37) 



(38) 



k<5. 



= -te l 



Using the relations ^ we write 



Jk(8 2 — 8i 



lk(* 3 -«0 



f 

— e 

t 



-3ik<5i 



-te 



i 



tk2 



-ik± 



l 



, //,•_ 



t 



-tc 



(k), 



(39) 



where we introduced the notations k\ = kai , hi = ka2 . Note that these variables 
change from — 7r to it. In the complex plane 0(k) delineates a circle centered at 
Co = 1 + e lk2 with a radius r = yT+ 2(<'/f) cosfc2 + (f /t) 2 when k\ is changed 
from — 7r to 7r for a fixed value /c2- In order to have a Dirac cone, the circle must cross 
the origin of coordinates. This is guaranteed for the values —3 < t'/t < 1. Thus, the 
Dirac cone is stable against additional interactions, diagonal hopping, when the last 
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condition is satisfied. The appearance of gapless Dirac fermions is not accidental to 
the honeycomb lattice, but is rather generic for a class of two-dimensional lattices 
that interpolate between square [t' — t) and 7T-flux (t' = — i) lattices. This indicates 
a certain robustness of the topological quantum number. 

3. QED2+1 description of graphene 
3.1. Noninteracting Dirac Hamiltonian 

Around two inequivalent K± points, where e(K±) = 0, the function (|11[) can 
be expa nded as 4>(K± + p) = ±Hvf(pi T V2), so that the Hamiltonian p0[) is 
linearizecP^ 



H ^ = Hl 1% [Tt(K + + P )H K+ (p)T (T (K + + p) 

~ J DC V Z7r ) 



(40) 



•P)] 



+ t;(k_ + p )w k _(p)t ct (k_ 

with the Hamiltonian densities for K± points 

H K + = hv F (ripi + T2P2), = hv F (-npi + t 2 P2)- 

The integration in Eq. (|40]) is done over the Dirac cone (DC) and the energy cutoff 
which preserves the number of states is 



(41) 




2.33i, 



(42) 



where Vt B = (2ir) 2 /S is BZ area. 

It is convenient to use a 4-component spinor ^(p) made from the two spinors 
for K± points. When the spinors T CT (K± + p) are combined into one, we exchange 
the sublattices in the spinor for K point @ so that 



**(p) 



^K+,<r(p) 
■0K_,<t(p) 



/a CT (K + 
ba(K+ - 
MK- 

Accordingly, the Hamiltonian (|40|) acquires the form 

K+A K+B 

( P X ~ iPy 

Px + iPy 



Ho = hv F 22 



DC 



dp 

d 2 p 



*1(p) 



V 



p 
p 

P)J 



K-B 







(43) 



\ 





-px + iPy 



*«r(p) 



/ 



= £ / 7^W^(p)^o(p)* ff (p), 

„ J DC ( 27r ) 



~Px ~ iPy 

W (P) = c^Pi + o?p 2 . (44) 



a This allows to use the same Hamiltonian for both K points, i.e. Hk ± = ihvp(ripi + T2P2) and 
obtain the standard representation of gamma matrices. 
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Here the 4x4 matrices a 1 ' 2 are the first two a matrices out of the three a 1 (i — 
1, 2, 3) matrices 

a = (a 1 , a 2 , a 3 ) =f 3 ® (ti,t 2 ,t 3 ) = (Z _M (45) 
that anticommute among themselves and with the matrix (3: 

/9 = n®l 2 =(°^ a ). (46) 

The irreducible representation of the Dirac algebra in (2 + 1) dimensions is given by 
2x2 matrices and there are two such irreducible representations differing in sign. 
They are both used in Eq. (|44|) , reflecting the fact that in addition to two degrees of 
freedom associated with the A and B sublattices, one should also take into account 
fermions at two distinct K± points. These degrees of freedom are described by 
t matrices that act in the valley (K±) subspace and as a result, we use a 4 x 4 
reducible (in 2 + 1 dimensions) representation of the Dirac matrices. 



3.2. Lagrangian of graphene 

To finish the mapping of the tight- binding model for graphene into QED2+1, we 
introduce 7 matrices and include an external magnetic field. We choose 

7° = A 7 = Pol = -if 2 « (ri,T 2 ,T 3 ) = f ° ~ T ^) (47) 

and introduce the Dirac conjugated spinor ^ , (T (p) = v E , j T (p)7°. The matrices 7" with 
v = 0, 1, 2, 3 satisfy the usual anticommutation relations 

{7^,7^-23^/4, <?^ = (1,-1, -1,-1), m," = 0,1,2,3, (48) 

but we remind that in QED2+1, we do not use 7 s in the Dirac representation of the 
Hamiltonian 

H * = E / 7^i(p)«?(p)^(p). 

~ J DC ( Z7T ) 

K (P) = 7°^o (p) - ftoF (7V1 + 1 2 P2 )■ (49) 

The representation (|47f is called the Weyl or chiral representation of 7 matrices 
and because these matrices are 4x4, one can construct the chiral 7 s matrix as 

„5 _ „-„ ; o„,i„,2„,3 _ Ma 



7 5 = nV7V = ^ _} 2 ] (50) 

which commutes with the a 1 and anticommutes with the 7". We note that the chiral 
representation Eq. (|47[) by means of the matrix 
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translates into another representation 

7" = fa«(r a ,ir 2 ,-ir I ) = (^ °),(^ ° ) , (52) 

which was used in Refs. |8j [T7J EH [55l Note that various researchers used different 
representations of the Dirac algebra (see e.g. Refs.lT|l6ll9ll56"|). so that the compar- 
ison of the U(4) symmetry breaking terms in terms of the underlying lattice and 
band structure of graphene is rather complicated. Our choice of the chiral represen- 
tation of 7 matrices in Eq. (|47|) is motivated by the fact that in this representation, 
the spinors with a definite chirality are eigenstates of j 5 which makes the physical 
interpretation of the chirality quantum number in graphene in Sec. 13.41 more trans- 
parent. Another advantage of this representation is that one can make a comparison 
with some recent works in Refs. 1571 1581 15T)1 1501 |6"T1 I61Z1 where the discrete symmetries 
of graphene were discussed. 

One can obtain from Eq. © the total electric current 

x>« =°L (w£ n(k) aH e T) T - (k) (53) 

where in the second line, we expanded the matrix near the K± points and intro- 
duced 4-component spinors ^'cr(p) and ^ , (T (p). One can recognize that the form of 
the electric current operator and its coupling to the vector potential A in Eq. © is 
standard for QED. Another famous condensed matter system, a d-wave supercon- 
ductor, is also described (see e.g. Refs. l63l l64l l65l l66l l67l [68)) in terms of QED 2 +i, 
but the coupling to the vector potential is different due to the presence of super- 
currents. 

Thus putting together the Hamiltonian, Eq. (|49p . and the interaction term 
A(n)j(n) (see Eq. (|6])), where current density j(n) corresponds to the current op- 
erator in Eq. (|53p . we arrive at the QED2+1 Lagrangian 

C = r ) [*7°(fift ~ Wo) + ihv Fl x D x + ihv Fl 2 D y ] r) (54) 

<T=±1 

with D a = d a + ie/hcA ai a = x,y written in the coordinate representation. In 
what follows we consider the setup when the external constant magnetic field B = 
V x A is applied perpendicular to the graphene plane along the positive z axis. The 
magnetic field also enters the Zeeman term included via \i a . When the field B is 
not perpendicular to the plane, the total field still enter the Zeeman term, while 
the orbital term would depend on its perpendicular com ponent . 

The theoretical explanation of the basic experiment J14 | 15 ] wri i cri proved that 
Dirac quasiparticles exist in graphene is grounded in the Dirac Lagrangian (|54[) . 
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Although the low-energy quasiparticle excitations described by (|54[) are nonintcr- 
acting, this Lagrangian captures the Dirac nature of the quasiparticles. M oreo ver, 
for magnetic fields below 10 T, there is no need to include the Zee man te rrrP^ and 
the orbital effect of the magnetic field is sufficient. Nevertheless, and acP^ 1 

measurements in strong fields B > 10 T already indicate that the interaction be- 
tween quasiparticles in graphene plays an important role. From a theoretical point of 
view, the effect of the long range Co ulomb in teraction had become a research topic 
long before the discovery of graphene (SEISED The Coulomb interaction Hamiltonian 
is 

Hc = ^fY.j rf 2 rrf 2 r'§,(r) 7 * CT (r)^-* CT ,(r')7 ^(r')- (55) 

a, a' 

The coupling constant is g = e 2 /ehvF = ac/evF, where a ~ 1/137 is the fine- 
structure constant. With air on one side of the graphene plane and SiC>2 on the other, 
the unscreened dielectric constant of the medium is estimated in Ref. [Ml to be e « 
1.6eo- The corresponding value of g « 1.37 is well above a simplifying assumption for 
theoretical analysis g <C 1 and this explains why the effect of Coulomb interaction 
in graphene is a hot topic of ongoing research (see e.g. Refs. [70l PfTl [72]) . 

Besides long range Coulomb interaction, one can also cons ider ot her lattice in- 
teractions like on-site repulsion and nearest neighbor 

repulsiorPlMI 

Hint = I E Tt(n)P + T ff (n)Tt,P + (n)T CT ,(n) (56) 
+ Y E Tt(n)P + T CT (n)T+,(n + 5 i )P-T CT ,(n + ^) 

n, Si ,o , ,a / 

with P± = (1 ±T3)/2. In the l ow-ene rgy continuum limit, this leads to several local 
four-fermion interaction terms E5|56| w j 1 j c j 1 m general break initial U(A) invariance 
and the generation of corresponding gaps (see, Sec. 13 . 6p : 



H mt =Y, U - ! d 2 r{^(v)Y^(v))\ r. ( = (P 7 3 , 7 5 ,7 3 7 5 )® ^0,^3). (57) 

i 

In summary, we would like to emphasize that the statement that the effective 
low-energy theory of graphene is massless QED 2+1 is based on three nontrivial facts: 

(i) Low-energy excitations in graphene are massless quasiparticles with linear dis- 
persion which have positive and negative branches, ±/juf|p|. It is often under- 
emphasized, that even the massive quasiparticles (see Sec. 13. 61 and Fig. 0(b) be- 
low) with the energy ±y/h 2 v F p 2 + A 2 are still governed by the massive QED2+1 
rather than Schrodinger theory. 

(ii) A qualitatively new feature of graphene is that the eigenfunctions of the low 
energy quasiparticle excitations obey the Dirac equation. As we have seen, the 
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spinor structure of the wave functions is a general consequence of t he hon eycomb 
lattice structure of graphene with two carbon atoms per unit cellPSS 
(iii) It is crucial that interaction of the quasiparticles in graphene with an external 
electromagnetic field be introduced using the minimal coupling prescription of 
quantum field theory. 



3.3. Continuum symmetries of QEDi+i model 

In the absence of Zeeman splitting, u a = u the effective Lagrangian ((54)) possesses 
global U(4) symmetry which is discussed for example in Appendices A and C in 
Refs. [9j [55l This symmetry operates in the valley-sublattice and spin spaces. It is 
useful to ignore, for a moment, spin space and begin by considering a global U(2) 
symmetry for the 4-component spinors in the valley-sublattice space. 



3.3.1. U(2) valley-sublattice symmetry 
One can easily check that matrices 



1 



1 



-I'Y 



-T2 ® T 3 , 



To 



1 



1 



;7 



-J3 <S> T , 



To 



1 



1 



fi <g> r 3 (58) 



2 ' 2 2 ' 2 ' 2 ' ' 2 

commute with the Hamiltonian ([44]) and satisfy the commutation relations of SU(2) 
algebra 

[Ti,Tj] = ie i]k Tk, (59) 

where 6y/- is the Levi-Civita sy mbol. Together with the identity matrix (T = I4/2), 
they lead to the U(2) syrnmetry^^ which acts in the space of valley and sublattice 
indices. We will call this effective symmetry of the low-energy approximation for the 
lattice graphene Hamiltonian as chiral symmetry because it contains the 7 s matrix 
and resembles chiral symmetries for massless particles in high energy physics (see 
a special Sec. 13.41 on chirality below). Because of chiral SU(2) symmetry, there is 
a conserving quantum number chirality which is characterized by eigenvalues of a 
diagonal generator T2 = 7 5 /2. Since there are two eigenvalues +1/2 and —1/2 of the 
diagonal generator T2, the four-component spinor a is the reducible representation 
of the SU (2) group. The conservation of chirality number plays an important role 
and leads to the absence of backscattering^] i n the presence of impurities that do 
not violate chiral symmetry. 



We note that the matrix 



S = 



/l \ 
10 




= s 1 =s~ 



10 

Vo -1/ 

brings the generators T,; to block-diagonal form 

l(r 2 0\ x l/r 3 

2U1J' ST2S -2{0r 3 



STxS- 



ST 3 S~ 



Ti 

n 



(60) 



(61) 
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In this reducible representation, the spinor \& a acquires the form 



/ a a {K++p) \ 
&.(K_+p) 
&.(K + +p) 

\-a a (K_+p)J 



(62) 



where the two upper and two lower components correspond to the two irreducible 
representations of SU (2) group. 



3.3.2. £7(4) spin-valley-sublattice symmetry and its breaking by the Zeeman 
term 

We now generalize the previous section and include rotations in spin space. The 16 
generators of the U(4) that operate in the spin and valley-sublattice space are 

^ 0/4 , fl®^, ^ 7 5 ? and |^ [7 3 l7 5], (63) 

where I4 is the 4x4 unit matrix and cr K , with k = 0, 1, 2, 3 being four Pauli matrices 
connected with spin degrees of freedom [<7o is the 2x2 unit matrix] . It is easy to 
see that when B = and there is no Zeeman splitting, the Lagrangian (|54|) and the 
interaction term (|55[) are invariant under global U (4) group generated by these 16 
generators. 

The simplest example of symmetry breaking is provided by the Zeeman term, 
^\a^ — ^ , 7 (T3^ / , where 173 acts on the spin indices a of the Dirac spinors. It 
explicitly breaks the E/(4) down to the U(2) a x U(2)t, with the generators 

°Y®h, ^f®il 3 , ^tt®7 5 , and ^®±[ 7 3 , 7 % (64) 

where k' = 0, 3. We will discuss other possible symmetry breaking terms in Sec. 13.61 
below. 



3.4. Chirality in QED 2 +i theory of graphene and its difference 
from the chirality in QEDs+i 

Because in the massless Dirac theory 7 s commutes with the Hamiltonian (f44|) (an- 
ticommutes with Hq given by (|49j) ). it introduces the conserving chirality quantum 
number which corresponds to the valley index. Indeed the spinors 

that describe the quasiparticle excitations at K± points, respectively, are the eigen- 
states of 7 5 : 

7 5 *k+=*k+, 7 5 *k_=-*k_- (66) 
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The labeling of K± points by eigenstates of the chiral operator 7 s is an advantage 
of the chiral representation (|47|) of the 7 matrices. This is not the case in the often 
used basis (|52|) with diagonal 7 . 

For massless particles in 3 + 1 dimensions, the chirality quantum number corre- 
sponds to the helicity which characterizes the projection of its spin on the direction 
of momentum. In the 2 + 1 dimensional case, the usual helicity concept for massless 
particles is meaningless and one may only talk about pseudohelicity. 

Let us now illustrate this by considering the Dirac equation which follows from 
the Lagrangian (|54[) in the simplest B = fi = case, also dropping the spin index. 
We consider its positive and negative energy solutions with a definite chirality 

*k ± M = e-^+^U^ ± (E,p), *k ± (t,r) = e l ^+^U^jE,p) (67) 

with E = Hvp\p\ which correspond to the electrons and holes from K± valleys, 
respectively. Substituting (|67p in the Dirac equation, we obtain that the spinors 

<(*.p) = (<f p) ). ^p)=(^_° E , p) ) («) 

satisfy Weyl equations, Ho(p)U e ' h (p) = ±EU e ' h (p) or for 2-component spinors, 

Hvf(tiPi + T 2 p2)V>K+ = i-EV'K+i 

- hv F {ripi +T 2 p 2 )ipli h = ±Eiplf_, (69) 

where the upper sign corresponds to the electrons and lower to the holes, respec- 
tively. We stress that in our case the vector p = (pi,p2) is in the graphenc plane. 
Formally, Eqs. look similar to the Dirac- Weyl equations that describe massless 
neutrinos)-^ 

ftcpa* = hc( P °' ° ) * = po*, (70) 

V -p°v 

but in the latter case, the space is 3D. This allows one to introduce, for a massless 
particle in 3 + 1 dimension, the helicity operator 

A =f s =(o:) ™ 

which commutes with the Dirac Hamiltonian hepa and characterizes the projection 
of particle spin on the direction of its momentum. Multiplying Eq. (|70[) by 7 5 
and taking into account that massless particles (antiparticles) have the dispersion 
Po = ±fic|p|, one obtains that 7 5 ^ = ±pS/|p|^. This illustrates that, for massless 
particles, the helicity coincides with the chirality, while for antiparticles it has the 
opposite sign to the chirality. 

Let us consider the solutions of Eq. (|69|) for the K + point 



1 / 1 \ ,. h >- - 1 




^kJ^p) = ^|I P,+ip, j, ^ + {E,p) = — y IpI j, (72) 
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and K_ point 

1 / 1 \ 1 / E2Ld£y\ 

V>K-(^P) = ^( P.+ip,j, ^AE,p) = —l IpI j. (73) 

The 4-component spinors made from the solutions l|72p. (|73p are mutually 
orthogonal eigenstates of the Hamiltonian H.q(p) and the operator 7 s : 

**(p)*2(p) = #**tb, ij = e,h, r, S = K+,K_. (74) 

In the 2 + 1 dimensional case one cannot make rotations around the direction of 
the quasiparticle momentum p lying in the 2D plane and the operator t = (ri,T2) 
does not have the physical meaning of the usual spin operator. Therefore the con- 
cept of the helicity for massless particles related to the Lorentz group and the real 

[7741 

space rotations is meaningless in this caseM- 1 ^ Furthermore, since there is only one 
generator of angular momentum r 3 , there is no non-abelian Lie algebra that can 
restrict its possible eigenvalues, giving rise to the possibility of exotic statistics (see 
Ref.[77|and references therein). The solutions ([72]) and ([73]) are double- valued under 
the rotations R{9) — exp(i8r 3 /2) with < 9 < 47r, because they describe spinor 
fieldsPl 

Formally one can consider the pseudohclicity operator 

.2D P S ^ ( T 



KZ =W S =loJ< T = ^ 2 ) (75) 

which commutes with the Hamiltonian and thus corresponds to a conserving quan- 
tum number. This operator, however, would correspond to an internal symmetry 
rather than to spatial symmetry as in the 3 + 1 dimensional case. Still one can relate 
chirality and pseudohelicity operators for massless quasiparticles, 7 s \& = ±A 2B ^. 

One can also see that the solutions in Eq. (j72|) and Eq. ([73]) which obey Eq. (l74l) 
differ only by the sign of momentum. Thus we arrive at the conclusior P^^^ 
that at the K + point, electronic excitations have energy ?ii\f|p| and Pt/IpIV'k = 
ipjc , while for holes the energy is —Hvf\p\, i-e. pT/\p\ip^ = — V'k • For the K_ 
point, these relations are inverted: for electrons, px/|p|-0| c _ = — rp^_ and for holes, 
Pt/IpIV'k. = V'k- ■ F° r 4-component spinors these conditions can be rewritten 
using the pseudohelicity operator: A 2D *| <:± = ±^^ ± and A 2D vI'^ ;± = zp^r^ They 
imply that at any given K± point, the direction of the momentum for electrons and 
holes with the same absolute value of the energy is opposite. This property of the 
massless electrons and holes is the consequence of equations of motion. However, to 
forbid backscattering of the quasiparticleJ^I one should also suppress the transfer 
of quasiparticles from one valley to another. This restriction is already associated 
with chirality conservation. For example, chirality is a good quantum number in a 
monolayer with electrostatic potential scattering. 

We were able to define the chiral matrix 7 s and the chirality quantum number, 
because we used the reducible 4x4 representation of the Dirac matrices. In 2 + 1 
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dimensions, there are two inequivalent irreducible 2x2 representations of the Dirac 
algebra!^ The corresponding 2x27 matrices are characterized by their signature®! 




One can check that 7" and matrices that constitute the 4x4 representation 
(I52p have the opposite signatures rj = +1 and 77 = — 1, respectively. It turns out 
that taking into account two inequivalent K.± points demands using two unitary 
inequivalent representations of 2 x 2 gamma matrices. As we saw in the chiral 
representation of 4 x 4 7 matrices, the K± points are distinguished by the chirality 
(valley) quantum number. The signature 77 of the corresponding 2x27 matrices is 

on 

sometimes known as "chirality" We note that by itself, this sign is not observable 
and one should introduce either the Dirac mass or/and the magnetic field and 
consider their relative signs (see Sec. 13.61 below). 

Finally, we note that in accordance with Eq. (I33p the matrix a 3 anticommutes 
with the Hamiltonian (|44l : {a 3 , Hq} = 0. The consequence of this is that the spec- 
trum of the Hamiltonian is symmetric with respect to the eigenvalue E — 0. The 
identity a 2 HQ(k)a 2 — —H (k) corresponds to Eq. ([34]) and guarantees the symme- 
try of the spectrum with respect to E = even in the presence of a finite Dirac 
mass , I'7 3 \E'. Also the matrix 7 anticommutes with the Hamiltonians (|4"4")l and AS)) : 
{7°,-£/o} = {7°!-ff(f} = 0- I n Refs. [4"6l I5U1 the property of anticommutativeness of 
Ho with a 3 matrix was called chiral symmetry, due to its analogy withanticom- 
mutativeness of 7 s with the Dirac Lagrangian in 3 + 1 dimensional casepHSSl Also 
sometimes the sign of the energy, e = Xvp\k.\, where A = 1 corresponds to the con- 
duction band and A = — 1 corresponds to the valence band is called the "chirality 
label" . We, however, keep the use of the word "chiral" for the characteristics which 
really involve the 7 s matrix. 

3.5. Discrete symmetries of QED2+1 model of graphene 

In the continuum description, the definitions of the discrete symmetry operations 
are not unique. In field theoretical studies of QED 2 +i, there is a certain convention 
on (see e.g. Refs.[l"2" ll79j) how to define these operations. It is based on the assertion 
that the parity transformation corresponds to inverting only one axis, say x axis: 
V(x,y) — ► (— x,y), because inverting both would be a rotation. Bearing in mind 
the condensed matter roots of the effective QED2+1 model, the discrete symmetry 
operations introduced here are in accord with the operations defined in Sec. 12.31 
for the tight-binding model. Knowledge of these symmetries allows us to properly 
classify possible symmetry breaking terms, avoiding the trap vividly described in 
Ref.il 

3.5.1. The spatial inversion V 

As discussed in Sec. 12.3.11 spatial inversion should invert both axes and exchange 
both A and B atoms and K± points. Thus applying the definition (|2"0|) and (|21[) for 




V = 



-TrfrVT 2 ]. 



(76) 
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the 4-component spinor ^ a (p) given by Eq. (|43]) , we defin^l 

*a(p) — '^(p)? -1 =P* ff (-p), P = f 1 ®T = 7°, P 2 = l. (77) 

Note that this operation of inversion of two spatial coordinates is not equivalent 
to a rotation through an angle tt in the plane which is given by ^(k x ,k y ) — -> 
exp(i7r(i7 1 7 2 /2)) , J , (— k x , — k y ) with ij 1 j 2 /2 being the generator of the correspond- 
ing rotation. 

The invariance (|22|) of Ho under V now follows from the condition 

PW?(p)P = «?(-p) (78) 

on the Dirac Hamiltonian density (j49| . Eq. (|78f generalizes the condition (f23|) on 
the 2x2 Hamiltonian density. It is easy to see that under this V transformation, 
the current ^uq^^ from Eq. (|53p transforms as follows 

^coT * -5" *o-o7 \l/, §tr 7 lv I / -^a-07 1 *, ^o^ 2 * Z -^a j 2 ^I, (79) 

where the unit matrix <7o acts on the spin indices a of the Dirac spinors. The 
interaction with fixed external magnetic field eA a ^UQ^ a '^ is invariant under spatial 
inversion. 



3.5.2. Time reversal T 



The time reversal operator, (|29|) interchanges K± points but not sublattices. Ac- 
cordingly when we define the action T on th e 4-co mponent spinors ^„ (p) , the part 
acting on the sublattice degree of freedom i J60 | 59 | r Jc| 



*(p) 



T = fl <g> 7} 



7 1 7 5 : 



i. 



(80) 



where <J2 acts on the spin indices of the spinor 4 , (J (p). Double T-transformation 
gives T 2 5'(p)T _2 = — 4 r (p). Hence for time-reversal symmetric systems with odd 
number of fermions, the degeneracy of levels cannot be less than 2 (Kramers' degen- 
eracy theorem). It is easy to see that K± points are indeed related by time- reversal 
symmetry 



*(p) = 



/o(K+ 
b(K + 
fe(K_ 

\a(K_ 



P)\ 
P 

P' 
p)J 



T 



T*(-p) = 



/o(K. 
6(K_ 
6(K 4 

Va(K, 



-P) 
-P) 

- P )7 



(81) 



where for simplicity, we ignored spin degrees of freedom. 



b The Vkj- components of the spinor (143 I t are exchanged by the inversion and because the sublat- 
tices are already exchanged by the definition of V'K. i we have to acting in the sublattice space 
instead of n used in Eq. (12 1 1 . This definition of V coincides with the corresponding operator in 
the second paper in Ref. 1591 

c Compare with the definition of V above which includes tq. 



February 1, _ 2008 11:55 WSPC/INSTRUCTION 

review-graphene'jversion 



FILE 



22 V.P. GUSYNIN, S.G. SHARAPOV, J. P. CARBOTTE 

In the time-coordinate representation, the transformation (I80|) and the corre- 
sponding transformations for the Dirac conjugated spinor are 

*(t, r) T*(t, r)T- x = icr 2 T*(-t, r), (82) 
*(t,r) - (r*(t,r)T- 1 ) t ( 7 )* = -*(-t,r)ia 2 T. 
The invariance (|30|) of Hq under T now follows from the condition 

Tng*(p)T = H£(-p). (83) 

on the Dirac Hamiltonian density (|49|) . Eq. (|83|) generalizes the condition l|3ip on 
the 2x2 Hamiltonian density. 

It is easy to see that under this T transformation, the current ^ero7 Al \I r from 
Eq. (|53|) transforms as follows 

vpo-cn 1 * ^ -^croT 1 *, ^ctoT 2 * -^cto7 2 *. (84) 



vpo^ * — > \1/(T 7 ^', 



When the external magnetic field is held fixed, the corresponding interaction term 
4'(7o7 Q ^ / now breaks the time-reversal symmetry. 



3.5.3. Charge conjugation C 

The charge conjugation C acquires nontrivial meaning when we consider the second 
quantized Dirac theory and introduce a transformation that exchanges particles and 
antiparticles (electrons and holes), leaving their spin and momentum unchanged. 
Accordingly we define the action of C on 4-component spinors ^(p): 

vMp) — ► C'MpJC- 1 = Cf^(p), C = 7 1 , (85) 

where T denotes the transpose. The matrix C satisfies the identitie^l 

C = -C- 1 = -C 1 = -C T . C^fC- 1 = -*f, n = 0,1,2. (86) 

It is easy to see that C exchanges A and B sublattices 



Mp) 



^a CT (K+ +p)^ 

& ct (k+ + P ; 
6 ct (k_ + p ; 

Va CT (K_+p)/ 
For the Dirac conjugated spinor we find 

(7VMP)) T 7° 



*«t(p) 



C*t(p) C "V 



/-6t(K+ + p)\ 
-4(K+ + p) 
4(K-+p) 
V 6t(K_ + p) / 



(-7 1 *.(P)) 5 



(87) 



It is easy to see that under C transformation, the current >I'ero7 Alv I' from Eq. (|53|) 
transforms as follows 

^0-07^*^-^07^*, n = 0,l,2. (89) 

One can also check that the Lagrangian (|54p is invariant under C if one simultane- 
ously replaces e — > — e. 



d Notice that our definition of C = 7 1 coincides with Ref. 1121 and that C(-y 3 ) T C 1 = 7 3 . 
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3.5.4. Difference between quasiparticles in graphene and massless neutrinos 

There is a widespread analogy between quasiparticles in graphene and massless neu- 
trinos which is based on the similarity of the Dirac-Weyl equations (|70|) and Eqs. (|69|) 
for quasiparticles in graphene. Our consideration of the discrete V, T and C sym- 
metries in Sees. I3.5TT| 13.5.21 and 13.5.31 shows that this analogy is not complete. The 
Dirac equation (|70|) corresponds to a pair of Weyl equations for 2-component left- 
handed, *l = (I/ 2 )! 1 - 7 5 )* and right-handed spinors ty R = (1/2)(1 + 7 5 )*. 
These two equations are related to each other by parity transformation. It is as- 
sumed in particle physics that in nature, there exist only left-handed neutrinos, so 
that it is sufficient to use one Weyl equation. However, thi s si ngle equation breaks 
C and V symmetries, but preserves CV and T symmetries !^ The Dirac quasipar- 
ticles in graphene with a definite chirality are also described by one of the Weyl 
equations for the 2-component spinor. The last equation, however, breaks V 
and T symmetries, but preserves VT and C symmetries. This conclusion about the 
difference of the discrete symmetries in QED3 + i and graphene could not be made 
by a formal comparison of Eqs. (|70[) and (f69|) . because one has to project the dis- 
crete symmetries of the lattice model from Sec. 12.31 into the continuum description 
of graphene. We also note that the fact that the K± points are related by time- 
reversal symmetry (see Eq. (181[) ) and this has an important implication for the 
Josephson effect in mesoscopic junctions consisting of a graphene layer contacted 
by two closely spaced superconducting electrodes! 2 -^ As Cooper pairs are made up 
of time reversed electron states, the two electrons in Cooper pairs that are injec ted 
from the superconducting electrodes into graphene go to opposite K± points!^ 
Finally, we note that the consideration of discrete symmetries allows one to make 
an analogy between the K± index and the circular polarization quantum number 
for photons. Indeed, photons with a definite circular polarization also break parity 
and time-reversal symmetries, but not charge conjugation. 

3.6. Dirac masses, their transformation under V , T , C and 
physical meaning 

One of the main reason why graphene is att racting the attention of many 
theoretician J 8 | 9 | 31 | 55 | 56 | 86 | 69 | 87 | ^^ | 89 | 90 | 91 | 92 | 93 J is that the electron-electron in- 
teractions and/or the magnetic field result in the breakdown of the f (4) symmetry 
and in new physics whose understanding demands going beyond the unconventional 
yet rather simple physics of noninteracting Dirac quasiparticles. 

For example, it is argued that the new zero filling factor state in the IQHE in 
graph ene, see n in a strong magne tic field B > 20 T is related to a spin polarized 
while v = ±1 states 26 28 are related to the lifting of sublattice or 
valley degeneracy. The spin polarized state is described by the order parameter 
(^ , 7°(T3\[') (or more general (fy^cr^)) which makes it similar to the Zeeman term, 
but in contrast, it originates from many-body interactions. 

Another channel of spontaneous breaking of the U(4) symmetry down to U(2) c x 
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U(2)d, is related to a possibility of the generation by the interactions of the Dirac 
mass Mo which enters the Lagrangian (|54|) in the following way 

C = *(t,r) «7° (ha d t - ia fj, + ia 3 ~(j, B B^ 

+ihv F <JQl 1 D x +ihv F <J -f 2 Dy - M D ] *(*,r). (90) 

The Dirac mass AfjJ has a general form Act k ®1 e Aa K T, where A is its absolute 
value, <7 K is one of the Pauli matrices (in what follows, we consider only <7q and 173, 
more general masses with cr were recently considered in Ref. 92) and T is one of the 
four matrices 

h, 7 3 , z 7 5 , 7 3 7 5 - (91) 

Under SU{2) group (as in Sec. 13.3.11 we ignore for a moment the spin degree of 
freedom) the gaps ^fPi>, \f'7 3 \E r , ^fij 5 ^ transform as a vector, while the gap ty^f 3 j 5 \f> 
is a scalar. Three gaps, ^/^,^j 3 ^, and ^ij 5 ^ break the SU(2) group down to £7(1) 
subgroup (with generators T3, T2 and T\ being unbroken, respectively). 

3.6.1. Physical meaning of the Dirac masses 

To understand the physical meaning of one of the masses, let us rewrite the Hamil- 
tonian Hi defined by Eq. (|24|) in terms of 4-component spinor ([43]) 

/ 7T -T I *t(p)[™+^®To + m_f 3 ®T3]* .(p), 
a J DC 

= / 7^2*(p)[™+^o7 + «i-ao7 3 ]*(p). (92) 
J DC ( 27r ) 

As was already mentioned after Eq. ([24]) . the term TO+CT07 changes the total carrier 
density of both A and B sublattices and can thus be absorbed in the term with the 
chemical potential /W07 . However, the term m_<7o7 3 breaks V defined by Eq. (|TT[) 
and is an example of one of the Dirac masses introduced above. The operator 

J d 2 x*a o7 3 * - Y,J DC -0j2 [4(K+ + PW(K + + p) + oj. (K_ + p)a ff (K_ + p) 

-bl (K+ + P )MK + + P) - &i(K_ + P )MK_ + p)] (93) 

determines the magnitude of the order parameter ($cto7 \?) (and hence the fermion 
gap) proportional to the electron density imbalance between the A and B sublattices 
and corresponds to the formation of a site-centered charge density wave (CDW) in 
the excitonic insulating ground stat which lifts sublattice degeneracy. As we have 
already mentioned in Sec. 12.3.11 it has been suggested^ that this kind of Dirac 
mass can be introduced in graphene by placing it on top an appropriate substrate 



e The energy gap Mjj is expressed via the corresponding Dirac mass m as Mo = mu|,. In what 
follows we ignore a difference between "Dirac gap" and "Dirac mass" and use the term "Dirac 
mass" for both of them. 
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which breaks the graphene sublattice symmetry between A and B and so generates 
an intrinsic Dirac mass for the fermions. 

The Dirac mass can also be generated dynamically. Historically the phenomenon 
of the electron- hole (fermion-antifermion) pairing in a magnetic field called magnetic 
catalysis was revealed in field theoryP^ Later, the experiments on graphite and their 
interpretation in terms of the Dirac fermioni ^] inspire d a theoretical condensed 
matter consideration of this kind of Dirac massP ^96 | 54 | Note ^hat in the notations 
of these papers, it corresponds to fyty. In the present conventions, the term Vf' 1 !' 
yields a chiral mixing between K± points: 

J d 2 zfa * = J2j Dc 4^2 [4(K+ + p)6<r(K_ + p) + bl (K+ + p)a CT (K_ + p) 
+fit (K_ + p)a a (K + + p) + 4(K_ + p)6 CT (K+ + p)] (94) 

This term preserves V , T and C symmetries. 

Magnetic catalysis is one of the candidates (see Ref. [31] for an overview of all 
scenarios) for a n explan ation of the new QHE states observed in graphene in high 
magnetic fields ! 26 1 27 1 2 ^ In particular, Refs. ESJ EH EQl EH are devoted to the fur- 
ther development of this scenario. The latest experimental r pose new 
questions that have to be addressed in the future theoretical work on spontaneous 
symmetry breaking and other competing models that attempt to explain the origin 
of v = 0, ±1 IQHE states. Here we will not go into the details of magnetic catalysis 
but instead consider all theoretically possible Dirac masses. Their physical meaning 
becomes clear when we write them in terms of the tight-binding model of graphene. 
Also we investigate their discrete symmetry properties. 

Similarly to the mass in Eq. (|94|) , the mass term with r = 7 5 also mixes K± 
points 

d 2 xVa h 5 y = -i Y,J [4(K+ + P)MK- + p) + (K+ + p)a CT (K_ + p) 
-fit (K_ + p)a CT (K + + p) - 4 (K_ + P )MK+ + P)] (95) 



/ 



and breaks V and C symmetries. The order parameters in Eqs. ([94)1 and ([95]) are 
related to a Kekule distortior|m w hich is ReA^^— ImA^ , i7 5 \E' and breaks in general 
V symmetry. 

Another interesting example is the mass term with F = j 3 j 5 for which one has 
J d 2 x*a o7 3 7 5 * = E/ 4^2 K(K+ + p)a CT (K + + p) - at (K_ + p)a ff (K_ + p) 
- fit ( K+ + p)MK + + p) + fit (K_ + P )fi CT (K_ + p)] . (96) 

In contrast to the gap (mass) in Eq. [M|) . the gap in Eq. ([M|) corresponds to a 
gap with the opposite sign at K_ point [cf. Eqs. ([M]) and ([9^]) ]. As pointed out in 
the second paper of Ref. [55] it is related to a model introduced by HaldaneSl as a 
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realization of the parity anomaly in 2 + 1 dimensional field theory (see also Refs. [3] 
[5l [83]) . Notice another definition of the spinor 

= (4(K+ + p), &+ (K+ + p), 4(K_ + p), 6t ( K _ + p)) (97) 

in Ref. [SSJ so that the gap in Eq. term becomes Vl^'Vo ® T3 ® T3\I/' (f and r 
matrices act in valley and sublattice spaces, respectively). In our conventions ([4"3"]l 
this corresponds to Vl'Vo ® fo ® = #<7o7 3 7 51 ^ gap- Comparing the masses in 
Eqs. (f9"5|) and ([TO]) , one finds that the second mass preserves V . On the other hand, 
considering that K± points are related by time-reversal T (see Eq. (|8ip). one can 
easily see that while the mass term A007 3 preserves T, the mass A7-<7o7 3 7 5 breaks 
it. This is also illustrated by the the corresponding energies of the LLL given by 
Eqs. (jlOip and (|106[) . respectively. In contrast, the mass ~ <737 3 7 5 which is related 
to the spin-orbit interactional does not break T. In fact, the gaps C737 3 7 5 and o$®I± 
are the only gaps that respect all the discrete symmetries. 

Yet another option for the Dirac mass which involves spin degrees of freedom 
was recently considered in Ref. [56] 

J d 2 x^a 3l 3 ^ = J2j 44i ° [4(K+ + PK(K+ + p) + 4(K_ + p)a a (K_ + p) 
- 4(K+ + p)b a (K + + p) - 4 (K_ + p)6„(K_ + p)] . (98) 

Since V transformation does not affect the spin variable, this mass breaks this 
symmetry as does the mass (j9"3"]l . However, in contrast to this latest case, the mass 
([98]) is also T breaking. 

All possible Dirac mass terms and their transformation properties under 

V, T, C symmetries are summarized in Table [1] For completeness we also included 
matrices M^, Mz and Ai^i , that correspond to the chemical potential, Zeeman 
term and generalized chemical potentials, respectively. For convenience of compar- 
ison with Ref. [58] we provide equivalents M' of these terms for $ ' and spinors 
([97]). 

Concerning the continuous E/(4) symmetry group, all Dirac gaps except <7o<g)7 3 7 5 
break it down to U(2) ® U(2) subgroup though the pattern of breaking depends on 
a concrete gap. For example, in case of the gap A^-f 3 ^ we have breakdown of U (4) 
to the U a {2) ® Ub{2) with the generators 

a a a a 

Y®/ 4 , y®r 2 . (99) 

Subgroups U a .b{2) act in the state spaces (K+,T), (K+,1) and (K_,f), (K_,l), 
not mixing these two spaces. 

An interesting example of U (4) symmetry breaking is provided by the gap 
\|r<7 3 ^/ 3 i|/. I n this case, we have 8 unbroken generators 

a 1 a 2 a 1 a 2 ct° a 3 <j° a 3 , , 

—®Ti, —®T U —®T 3 , —®T 3 , —®T 2 , —®T 2 , —®h, — ®/ 4 ,(100) 
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Table 1. Transformation properties of the various Dirac mass terms 
^Mrj* = *tM* = *'tM'*' with M D = ero,3 ® T under 7>, T and C trans- 
formations. The matrices Mp, Mz and A^ui,«Mu2 correspond to the chemical 
potential, Zeeman term and generalized chemical potentials, respectively. 
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which form another U a ' (2) ® Uy (2) group acting in the state spaces (K + , |), (K_, J.) 
and (K+, 1), (K_, |), respectively. While this gap is "P-odd, T-odd it is invariant 
under combined 'PT-inversion. 

Finishing our classification of the Dirac masses, we mention that the 
QED2+1 desc ription of d-wave superconductivity also involves various Dirac 
masse j£- 4 | 65 | 66 | 67 168| w ] aen the possible opening of a secondary gap is considered 
(for example, a charge or spin density wave or a second superconducting order 
parameter with different symmetry). 



3.6.2. Dirac Landau levels 

In an external magnetic field B Dirac Landau levels are formed. In the presence of 
a finite gap amplitude A ^ and ignoring the Zeeman splitting, for the matrices 
T = I4, 7 3 , 7' 5 the energies of Landau levels are 

r±Asgn(eB), n = 0, 

" \sgn(n)y/A 2 + 2\n\hv 2 F \eB\/c, n = ±l,±2,.... 1 ' 

To understand how the spectrum (|101[) emerges, it is instructive to write the Dirac 
equation for one of these masses, considering separately each K point. For example, 
for the A(T 7 3 mass, one obtains 

[irohdt + ihvpT\D x + ihv F T 2 D y - t 3 A] vp K+ (t,r) = 0, (102) 
[iT hd t — ihvFT\D x — ihvFTiDy + t 3 A] ?/>k_ (t, r) = 0. (103) 
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Fig. 3. (Colour online) Band structure of graphene. Electrons are shown in red and holes in blue, 
(a) The low-energy linear-dispersion -E(k) near the Dirac K+ and K_ points for B = 0. (b) A 
possible modification of the quasiparticle spectrum by the finite gap (Dirac mass) A. The chemical 
potential (indicated by horizontal line) /i is shifted from zero by the gate voltage, (c) Landau levels 
E n in the Dirac theory of graphene. Spin degree of freedom is ignored. For a given direction of the 
magnetic field B applied perpendicular to graphene's plane, the lowest (n = 0) Landau level has 
the energy Eq = —A at K+ and Eq = A at K_. 



Multiplying the left side of Eq. ([TU2|| by t 3 and the left side of Eq. |(TD5|) by -r 3 , 
one obtains 

[ifhdt + ihvpl 1 D x + ihv F fDy - A] Vk+ (i, r) = 0, (104) 

[ifhdt + ihv F ^ x D x + ihv F j 2 D y - A] ip K _ (t, r) = 0, (105) 

where in Eq. (|104[) . 2x2 gamma matrices are given by Eq. (|52[> and in Eq. (|105[) 
7° = — 7°, 7 ' = 7 ' and have the opposite signature rj. For reference, the solution 
of Eqs. (|104|) and (|105[) in an external magnetic field is written down in Appendix A. 
The energies of the Landau levels with n and — n for |n| > 1 are symmetric (see 
Eqs. (|101[) and (|106p ). The LLL is, however, asymmetric and anomalous, because 
for Eq. (|104p its energy is — Asgn(eS), while for Eq. (|105[) . its energy is Asgn(eS) 
(see Fig. 0(c) where the case eB > is shown). It turns out that the sign of the LLL 
energy is defined by the relative sign of the signature r\ of the 2x2 Dirac matrices 
(|75|) and the Dirac mass A^Hl This implies that because the sign before the mass 
A in Eqs. ()104|1 and () 1 05(1 is the same, the sign of the LLL energy depends on the 
sign of the signature of 7" and 7^ matrices. Thus altogether the combined energy 
spectrum of Eqs. (|102p and (|103|) is symmetric and given by Eq. (jlOip . This reflects 
the fact that the Dirac mass Acto7 3 does not break the time-reversal symmetry. 

In the limit A — > 0, the Landau levels ±Asgn(ei?) corresponding to n = 
Landau levels merge together to form a single level. However, for fj, = and ignoring 
Zeeman splitting, it remains half- filled. This property of the LLL equally shared by 
particles and antiparticles is at the heart of the unconventional QHE in grapheneP^ 

For F = 7 3 7 5 , the spectrum is asymmetric 

f-Asgn(eB), n = 0, 

" \ sgn(n)-v/A 2 + 2\n\hv F \eB\/c, n = ±l,±2,... 1 ' 

Indeed, the T breaking Arco7 3 7 5 mass results in two equations 

[iT hd t + ihv F T\D x + ihv F T 2 D y - t 3 A t ] V>k + (*> r) = 0, 

[irohdt — ihv F T\D x — ihv F TiD y — T3A7-] V'k. (t, r) = (107) 
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which have the same sign before A7-. Accordingly, because the corresponding 2x2 
gamma matrices are the same as for Eqs. (|104j) and (|105[) . one can readily see that 
both equations in (|107p lead to the same sign of the LLL energy, viz. — At sgn(eS). 
Accordingly, the LLL remains doubly degenerate. Therefore, the combined energy 
spectrum (|106j) is obviously asymmetric. Thus, the time-reversal breaking character 
of the Dirac mass Atoo7 3 7 5 is revealed by an external magnetic field. At temper- 
atures well below At this leads to the QHE even in zero magnetic field.^ Indeed, 
since the transverse Hall conductivity a xy is odd under time reversal, a nonzero a xy 
can occur if time-reversal invariance is broken. 

In the limit A3 — ► 0, the n = Landau level — Asgn(e_B) in Eq. (|106|) is either 
electron- like or hole-like depending on the relative sign of A and eB. Thus, the 
generic zero gap limit is indefinite in terms of a clear particle or hole character 
of quasiparticles occupying the LLL. Nevertheless, from a physical point of view, 
the property of the LLL to be equally shared by particles and antiparticles can be 
argued by the presence of the Zeeman term. 

The degeneracy of Landau levels with n = ±1, ±2, ... in Eqs. (| 1 1 [) and (|106p is 
\eB\/(irh) (per unit area and per spin), while the degeneracy of the lowest Landau 
levels (LLL) n = for Eq. (fIOT|) is \eB\/{2irh) and for Eq. (TiTMj) is \eB\/{irh). 

It is worth emphasizing that it is the spectacular "relativistic" energy Landau 
scale L(B) = ^J\eB\hv 2 F /c which defines the distance between Landau levels 

AE = E Y (A = 0) - E (A = 0) = ^/2L 2 (B) w 424v/£[T]K. (108) 

For instance, it corresponds to AE w 2800 K at B — 45 T which makes it possible 
for the QHE in graphene to be observed at room temperaturePSl 

Including the real spin degree of freedom, we have the following situation for 
Landau levels in a magnetic field with the "P-odd, T-even gap For levels 

n > 1 all Landau levels are 4-fold degenerate (states |K+,f), |K + ,j,), |K_,t), 
|K_, J.)). On the other hand, for the n — Landau levels, as shown in Fig.|4](a), we 
have 2-fold degeneracy: states |K+, f), |K+, j) and |K_, f), |K_, J.) with energies 
=FA, respectively. For V-odd, T-odd gap A^a^ 3 ^ (see Fig. @] (b)) the situation 
with levels n > 1 is the same, but two Landau levels with n = with energies 
=FA now have different 2-fold degeneracies: |K+,T), |K_,1) and |K_,f), |K + ,1) 
corresponding to a symmetry breaking pattern described by Eq. (|99p . In both cases, 
switching on the Zeeman interaction leads to further splitting of the level n = 
resulting in four different levels, thus completely lifting the 4-fold degeneracy of this 
level in noninteracting theory. For both gaps considered we expect QHE with the 
filling factors v = 0, ±1, ±2fc where k is a positive integer. 

Let us now consider the P-even, T-odd gap A^'7 3 7 5 \E' which does not break 
the 17(4) symmetry. All Landau levels have 4- fold degeneracy. More precisely, for 
the levels with n — 0, the states of either spin from K + - valley residing on the 
sublattice B are degenerate with the states of both spins from K + - valley residing 



We are grateful to D.V. Khveshchenko for a discussion on Dirac masses considered in this section. 
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Fig. 4. (Colour online) Landau levels including spin degree of freedom, but taking zero value of 
Zeeman splitting. Electrons are shown in red and holes in blue, (a) for <tq <8>7 3 gap. (b) for CT3 ® 7 3 
gap. 



on the sublattice A. Applying the Zeeman term, one now obtains only 2 different 
(doubly degenerate) levels. 

Finally, for "P-even, T-even gap /S.^>a^ 3r y b ^ one has two pairs of states for the 
n = Landau levels: |K + , f), |K_, f) and |K + , J.), |K_, J.), which remain degenerate 
even when the Zeeman interaction is turned on. Hence, for the last two gaps we 
expect QHE with the fillings factors v — 2k with k — 0, ±1, ±2, . . .. 



4. Electromagnetic response 
4.1. Electrical conductivity 

The frequency-dependent electrical conductivity tensor a a p(VL) is calculated using 
the Kubo formula 

K aP (Q + iO) (r a g) n&(ft + fO) 
-.(fi + zO) ' ^(« + »0)=— + ^ , (109) 

where the retarded correlation function for currents is given by 

/oo 
dte int n^(t), U^(t) = -ie(t)TT(p[J a (t),J p (Q)]), (110) 
-00 

V is the volume (area) of the system, p — ex-p(—/3Ho)/Z is the density matrix of the 
grand canonical ensemble, (3 = 1/T is the inverse temperature, Z = Trexp(— (3Hq) 
is the partition function, and J a are the total paramagnetic current operators with 

J a (t) = e im ' n J a {Q)e- im ' n , J a (t) = 52£&n), (HI) 

n 

expressed via the paramagnetic current density (jSJ) . The diamagnetic or stress tensor 
(r a( g) in the Kubo formula (| 109[) is a thermal average of the diamagnetic part ([9]) 



r a (3) 



(E^( n ))- ( 112 ) 



In many cases, the effective low-energy QED 2 +i description provides a very good 
starting point for the investigation of the various transport properties of graphene 
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making it possible to use the powerful field theoretical methods that allow us to 
obtain rather simple analytical expressions. In this case, the conductivity tensor in 
Eq. (|109j) is calculated for the model Lagrangian ([90]) rather than the full tight- 
binding model in Eq. (|5j). 



4.2. Zero field AC conductivity in the continuum Dirac model 



One of the interesting results for the zero magnetic field case is that in the 
high-fre quency limit, the interband contribution to longitudinal conductivity is 
constant ffllOOflOT] 



2 

7TC 



Rea xx (fl) ~ — , 0>|/x|,T. (113) 

It is remarkable that the conductivity is universal, independent of the band structure 
parameters t an d vp . Moreover, Eq. (|1 13|) is even valid for a finite wave-vector k, 



kvp <C T <C fiffSOl The flatness of the conductivity results from the interband 
transitions between electron and hole bands at K± points. Of course, this result 
needs modification if the photon energy becomes comparable to the band edge 
energy, when the assumption of linear momentum dispersion relation ceases to be 
valid. We will return to this important question in Sec. 14.41 

A more general expression^ for the conductivity cr xx (Cl, T) valid in the limit 
of small impurity scattering rate T(u>) and neglecting the real part of the impurity 
self-energy reads 



a xx {ft,T) = — / dto 



[rip (w) — nir(w')] 7r 



n 2 + 4T 2 (uj) (uj + oj') 2 +4r 2 (w)_ 
x (M + \lu'\)(uj 2 +uj' 2 ), uo'=uj + n. (114) 

The first term in square brackets of Eq. (|114[) describes the intraband transitions 
and the second term, which reduces to Eq. (|113|) .describes the interband transitions. 
An essential feature of Eq. (|114p is that we kept the energy dependence of T(u>). In 
deriving this equation, we have assumed the small Q limit, so that r(a/) ~ r(w) 
and for Q <C T, the difference [rip (w) — np(u>')]/Cl can be replaced by the derivative 
—dnp(uj)/diu. The intraband term of Eq. (|1 14|) in this case results in the following 
expression for Drude conductivity 

^(n,T) = a 00 JdJ^^) >^\ , (115) 



duj j n 2 + AT 2 {u 

— oo 

with coo = e 2 /(ir 2 h). As discussed in Ref. [551 the frequency dependence of T(uj) 
could be probed by measuring a xx (fi) at the different values of the gate voltage V g . 
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In fact, if it is assumed that T(lo) is linear in u>, as can be expected in the 
Born approximation for weak scattering when the chemical potential fi = 0, i.e. 
r(w) = 700 + ct\u>\ with a small value of 700, then one can show 

7T tt~ 



la 



1 8a T 



700 <!i«T. (116) 



Note the linear dependence on microwave frequency with slope inversely propor- 
tional to the temperature. A similar formula has been derived for the microwave con- 
ductivity in a d-wave superconductorp^H a system which, as mentioned in Sec. 13.21 
can also be described by QED2+1. It was used to explain the cusp like behavior seen 
in pure samples of ortho II YBCO6.5 in which every second chain is complete and the 
others entirely missing. Eq. (|115p offers the possibility of observing directly through 
microwave experiments the transport quasiparticle scattering rates in graphene and 
establishing whether the scattering is weak and the Born approximation applies, 
or it is strong (unitary limit) and a better model for T(w) is constant Tq within 
the energy range of importance in microwave experiments. In this case, Eq. (|115p 
reduces to a Drude profile even for /i = 0. 

Another prediction of theory is that for > T, Eq. (|1 16(1 is replaced by 

a xx (n,T)^a 00 2^\ ^ 0a+a ^ l) (117) 
" + 4(700 + a\fjt\y 

which has the Drude form and shows the remarkable property that its width can be 
increased continuously by increasing the gate voltage, although the impurity content 
is left unchanged. This arises due to the relation [i oc sgn V g y/\V g \ (see Sec. 14.4.21 
below), because of the increase in residual scattering which is proportional to the 
final state density of states which varies as |w|. 

In the context of massive Dirac quasiparticles, we note that for T = and in 
the limit of zero scattering T — > 0, one can obtain a remarkably simple form for the 
conductivity, viz. 

a „ (0 ) . y <(n) ^-Ay-A') (118) 



2h n 2 V 2 



Here, the theta function 6 (|fi|/2 — max(|/i|, A)) cuts off the low f2 part of the inter- 
band contribution at 2|/i| or 2A whichever is the largest. The expression in Eq. (|118p 
was obtained considering as an example the excitonic gap in Eq. (|93[) . The factor 
(£1 2 + 4A 2 )/£! 2 in the second term of Eq. (|118jl modifies the interband contribution 
from its constant value of the massless Dirac case. At large ft, this modulating factor 
goes to 1, but for = 2 A it is equal to 2. This shows that the existence of a finite 
A does have a clear signature in the AC conductivity. So far only ARPES measure- 
ments show some evide nce f or a finite value of the gap A at B = in graphene 
epitaxially grown on SicP^ We will consider the situation for finite magnetic fields 
in the next Section. 
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4.3. Magneto -optical conductivity 

Motivate d b y recent experimental advance s in infrared spectroscopy of one layer 
graphen eP^ a few layer epitaxial 

gra p h itpl 

and in highly oriented pyrolytic 
graphite ] * we discuss some results obtained in Refs. 1106} l99l for magneto-optical 
conductivity of graphene here. Starting from the Lagrangian (|90j) . one can obtain a 
simple representation for the complex diagonal conductivity 

a xx m = -^M (119) 

°° ( / A 2 \ 

x l2 { [ X ~ M n M n+1 ) ^ UF ^ " n ^ M «+i)] + [n F (-M n+1 ) - n F (-M n )}) 



M n+1 - Mr, 



1 1 



M n - u n+1 + n + *(r„ + r n+1 ) M n - M n+1 -n- i(T n + r n+1 ) 

A 2 \ 

{[n F (-M n ) - n F {M n+1 )} + [n F (-M n+1 ) - n F (M n )]) 



M n M n+1 
1 



M n+1 + M n 



1 



M n + M n +i + n + i(r„ + r„+i) M n + M n+1 -n- i(T n + r„+i) 



where M n — _\/A 2 + 2nvp\eB\/c are the absolute values of the energies of the 
Landau levels|3 The Zeeman splitting was not taken into account, because it can 
only shift the positions of the lines if the value of the spitting depends on Landau 
level index n. The Lorentzian representation (|119p is derived in the Appendix of 
Ref. [99] using the assumption that the Landau level half-width r n depends only 
on the Landau level index n and is independent of the energy ui. In general, the 
scattering rate T n (uj) is expressed via the retarded fermion self-energy, T n (uj) — 
-ImE^'(w), which depends on the energy, temperature, field and the Landau levels 
index n. This self-energy, which in general has a real part also, has to be determined 
self-consist ently fr om the Schwinger-Dyson equation. This equation can be solved 
analytically^!^ and numerically, as was done in Refs. [TBI HB When the frequency 
dependence of ImS^(w) is retained, the ReE^(w) is nonzero, which results in the 
shift of Landau level positions. A n analy tical expression for the optical conductivity 
can also be obtained for this case! ' ' Similarly to Eq. (|118p . the expression (|119[) 
was obtained by in cluding the excitonic gap ([M]) . 

Jiang et a/M^ have recently performed optical transmission experiment in an 
external field on a single layer graphene. Earlier measurements were on several layer 
systems, namely ultra thi n ep itaxial graphitefESI Another recent work reported a 
magneto-reflectance stud y~^5 | f highly oriented pyrolytic graphite in a magnetic 



? The levels with the energies E n = —M n are explicitly included in Eq. (1 11 9 P and we set H = 1. 
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Fig. 5. (Colour online) Real part of the longitudinal conductivity, Hea xx (£l) in units of e 2 /h 
vs frequency Q in cm -1 for field B = 18 T, temperature T = 5K. Long dashed, the chemical 
potential fi = 1310K and the gap A = OK, dash-dotted fi = 2930K and A = OK, solid fi = 1310K 
and A = 250 K, short dashed fi = 2930K and A = 700 K. For the left frame the scattering rate 
r = 60K and for the right frame scattering rate T = 190K. 

field up to 18 T. In this experiment the Landau level energies are found to be linear 
in B, while in graphene as well as in the epitaxial graphite samples they were found 
to go like \f~B instead, which is consistent with Dirac quasiparticles. While this is 
expected in a single sheet graphene, the situation is not as clear for the multilayer 
case. Consequently the data in this case were analyzed in some detail in Ref. 11061 
and the discrepancies with our Eq. (|119|) were noted. In particular, it is hard to 
understand the observed ratio of the intensity of the first to the second interband 
li ne. W e do not repeat details here, rather we comment on the newest data of Jiang et 
a/P21. In particular, these authors find that the distance between the Landau levels 
with the energies —M\ and M2 is somewhat bigger than predicted on the basis 
of the independent particle model such as Eq. (|119p . They exam ine the possible 
effect of Coulomb interaction as described by A. Iyengar at a/ff^on Landau levels 
and tentatively conclude that this could possibly explain the data. On the other 
hand, they observe that the model with a gap for the v = plateau included as in 
Eq. (|119p is not compatible with their observations, at least in the present form. 
It predicts a reduction rather than the observed small increase in spacing between 
first and second interband optical line. At this early stage, however, it is hard to 
be definitive about correlation effects. In Fig. [5] we show results for the real part 
of the longitudinal conductivity Rea xx (fl) for B = 18 T, T = 5K and two possible 
Landau level half- width T = 60 K (roughly 4 times larger than that observed in dc 
measurements for v — 4 plateaux in Ref. [2"6"|) and T = 190 K, which corresponds 
to the observed (in Ref. [2"9"| 20 fs half- widths of the transmission peaks. We show 
results for two values of chemical potential p = 1310 K and n — 2930 K. The first 
corresponds to the line Ti in Ref. [551 measured at the filling factor v = 2 with 
the carrier imbalance p = 9.4 x 10 11 cm~ 2 , and falls between Mq and Mi Landau 
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level energies in the notation of Eq. (|119|) . The second corresponds to the line T 2 in 
Ref. [29l v = 10 with the carrier imbalance p = 4.7 x 10 12 cm -2 and falls between 
M 2 and M 3 Landau level energies^] The long dashed (red) curve includes no gap A. 
It displays an interband absorption peak at = 1248 cm -1 corresponding to the 
transition from n = to n — 1 levels, with a second at fi = 3013 cm -1 (from —Mi 
to M 2 and -M 2 to Mi), a third at fl = 3927 cm -1 (from -M 2 to M 3 and -M 3 
to M2) etc. The solid (blue) line is for the same case, but now a gap A = 250 K is 
introduced to illustrate its effect on the absorption line. We see no change in the 
higher energy peaks but the peak at il = 1248 cm -1 is split into two, one shifted to 
lower energies and the other to higher energies. Each one is broadened by the peak 
half- width 2r = 120 K. While our parameters were chosen with the experiments of 
Jiang et al.^ in mind, no splitting of the first interband line is observed, which is 
consistent with the fact that their T is considerably larger. If we increase T instead 
to their value 190 K, the results are shown in the right hand frame, which differs 
from the left hand counterpart only through the value of F. We see broader peaks 
associated with the various optically induced allowed transitions between Landau 
levels. Also, the broadening of the first interband line in the solid (blue) curve is 
now sufficient, so that a single peak is seen. However, it does not have the single 
Lorentzian lineshape of the dashed (red) line. There is a slight depression at the 
position of the original peak, a flat top and slight overall broadening of the line. 
There is no visible increase in the distance between the first and second absorption 
peak. 

We have included in Fig. [5j for comparison, results for another value of chemical 
potential fi = 2930 K. This corresponds to the case when \i lies between the second 
and third Landau levels. As described in Ref. 11061 in this instance the first interband 
line has faded into the background, as has the second, and the third has halved its 
original intensity (i.e. compared with the long dashed (red) curve). This is seen in 
the dash-dotted (black) curve which has A — 0. Note that a new line appears at 
lower energies Q — 397 cm -1 which corresponds to an intraband transition between 
the second and third Landau levels. To observe even a small change in this curve 
in the region above ~ 3500 cm -1 , we need to include a large gap A = 700 K as seen 
in the dotted (green) curve. 

One of the important predictions of Eq. (|119|) which has not been addressed 
yet in any of the reports on optical experiments done so far, is the behavior of the 
absorption lines as the chemical potential is changed. As we saw in Fig. [5] for /i 
falling in the energy range between n = and n — 1 Landau levels, there is a line at 
£1 = Mi . This line is anomalous in the sense that it disappears when fi falls beyond 
this range. If /1 falls between n — 1 and n — 2, there will still be an interband line 
at Mi + M2, but it will have only half the intensity it had in the previous case. The 



h The actual measurements in Ref. 1291 were made for holes with v = —2 and —10, but from 
theoretical consideration the particle-hole symmetry is obeyed because the conductivity ||119| I is 
an even function of fj.. 
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other lines remain. If we further increase u to fall between M 2 and M 3 , this would 
correspond to the dash-dotted (black) curve in Fig. [5] discussed above, i.e. the first 
and second interband lines have disappeared and the third halved its intensity. 

It is hoped that future experiments will verify this predicted pattern of 
change with incre ase of /i and in particular verify that the first interband line 
is anomalous ESI it appears at full intensity or not at all, while all other lines 
first fall to half intensity before disappearing entirely. Its optical spectral weight is 
transferred to an intraband line at M n+ \ — M n . 

We note that the present consideration is based on a perfect symmetry between 
particles and antiparticles, which is built into the Dirac formalism. Accordingly, as 
we saw above, two different tra nsit ions can contribute equally to the same optical 
line. A most recent experiment!^ shows, however, clear evidence for the break- 
ing of particle-antiparticle symmetry in the graphene system at the level of 2.5%, 
approximately five times larger than expected. 

Coming b ack to the role of correlation effects, we mention another recent paper 
by Jiang et alW^ which also points to possible limitations of an independent particle 
theory. Based on magnetoresistance measurements in fields up to 45 T it is concluded 
that the observed activation gap is much smaller than expected on the basis of the 
known spacing between Landau levels, which is set by the value of the Fermi velocity 
vp. These results point to a need to include in the theory, additional effects which 
go beyond the independent particle model. In this work, we have classified and 
discussed several possibilities related to different Dirac masses and, as an example, 
considered the effect of an excitonic gap on the optical conductivity. We hope that 
this can help when pursuing future developments. 

4.4. Sum rules for the optical and Hall conductivity 

It is instructive to consider the case when QED2+1 approximation is not sufficient 
and one should consider a full tight-binding model . This example is given by the 
optical conductivity sum rulesM^The real part of the frequency dependent optical 
conductivity er xx (f2) is its absorptive part and its spectral weight distribution as a 
function of energy (M7) is encoded with information on the nature of the possible 
electronic transitions resulting from the absorption of a photon. Even though the 
relationship of the conductivity to the electronic structure and transport lifetimes is 
not straightforward, much valuable information can be obtained from such data. In 
particular, the /-sum rule on the real part of cr xx (Q) states that (see e.g. Refs. l40l 



where (t xx ) is the thermal average of the diamagnetic term (|112[) . The optical 
conductivity sum rule is a consequence of gauge invariance and causality. Gauge 
invariance dictates the way that the vector potential enters Eq. ([5]) and determines 



[TTOl ITT21 ITTT1 for a review) 



OO 




(120) 



— OO 
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the diamagnetic and paramagnetic terms in the expansion (J6j) as well as the form 
of Kubo formula (|109[) . Causality implies that the conductivity, Eq. (|109[) satis- 
fies the Kramers-Kronig relation. Note that defining the plasma frequency up via 
ujp/(4ir)5 a p = (r a i3)/V, we can rewrite the RHS of the sum rule (|120ll in terms of 

Up. 

For a finite single tight-binding band, one can explicitly calculate the ther- 
mal average (11121) and obtain the following representation of the optical sum 
rul j40lll0lll2lllll 

dORea^fi) = _£n k>CT ^, (121) 

k , (7 x 

with Qm being a cutoff energy on the band of interest and only the contribution 
to Re a xx (£l) of this particular band is to be included in the integral. Here is 
the electronic dispersion, k is the wave vector in the Brillouin zone, and nk, CT is 
the probability of occupation of the state |k, a). For tight-binding dispersion with 
nearest neighbor hopping on a square lattice, it is easy to show that the right-hand 
side (RHS) of Eq. (|121[) reduces to e 2 /h 2 multiplied by minus one-half of the kinetic 
energy, Wk.e. per atom. 

When a constant external magnetic field is applied to a metallic system in the 
z-direction, the optical conductivity acquires a transverse component cr xy (Q) in ad- 
dition to the longitudinal component cr xx (p,). This quantity gives additional infor- 
mation on the electronic properties modified by the magnetic field. In this casep^ 
there is a new sum rule on the optical Hall angle 9h(Q)- If we define 

t H (n) = ten0*(n) = ^M, (122) 



then 



2 f°° 

- / dnRet H (n) = uj Hi (123) 
* Jo 



where u>h is the Hall frequency which corresponds to the cyclotron frequency lu c 
eB/mc for free electrons. 



4.4.1. Limitations of the Dirac approximation 

The fact mentioned above that the effective low-energy Dirac description of 
graphene is insufficient for the derivation of the sum rules can be easily under- 
st ood from two e xamples. As we have seen the RHS of Eq. (|120p is normally equal 
t(J 40 | 110 | lll | 112 fy k e thermal average of the diamagnetic term (|112[) that is defined 
as the second derivative of the Hamiltonian [see Eq. © above] with respect to the 
vector potential. This term is zero if one tries to obtain it from the approximated 
Dirac Lagrangian (f54"|) . On the other hand, it follows from the same Dirac approxi- 
mation t h at , in t he h igh-f req uency limit, the interband contribution to conductivity 
is constant! ^^^ ^^ given by Eq. (|113[) . The frequency fl is unbounded in the 



February 1, _ 2008 11:55 WSPC/INSTRUCTION 

review-graphene'jversion 



FILE 



38 V.P. GUSYNIN, S.G. SHARAPOV, J. P. CARBOTTE 

Dirac approximation, although physically f2 should be well below the band edge. 
This example indicates that when considering sum rules, one should go beyond the 
Dirac approximation. 

The second example is that the cyclotron frequency uj c for the Dirac 
quasiparticleJ^ is defined in a different way, lo c = eBv F /(c\fj,\), where vp is the 
Ferm i velocity . This definition follows from the fact that a fictitious "relativistic" 
mass SH l5 | 17 | m ^ _ plays the role of the cyclotr on mass in the temper- 

ature factor of the Lifshits-Kosevich formula for graphene 1 1 14 | 1 1 5 ] xhis cyclotron 
frequency diverges as \i — > 0, also posing the question as what one should use as a 
Hall frequency on the RHS of Eq. (fT23| . 

4.4.2. Diagonal conductivity sum rule 

These problems were resolved in Ref. I109[ where considering a tight-binding model 
|5]) we obtained the RHS of the sum rule (|120|) to be equal to 



2< J Bz ^^^))- n ^e{m 



V h 2 



d 2 _ ( d<p(k) \ 

Ok 2 ' \ dka J 



e(k). (124) 



Here njr(w) = l/[exp((w — y)/T) + 1] is the Fermi distribution, the phase (fi(k) is 
defined below Eq. (|10[) and the energy e(k) is given by Eq. (|12p . The momentum 
integration in Eq. (| 1 24[) is over the entire BZ and the thermal factors rip (e(k)) and 
np{— e(k)) refer to the upper and lower Dirac cones (see Fig. [5]), respectively. We 
stress that a simple generalization of Eq. (|121[) for a two band case would miss the 
term with the derivative of the phase, (dip(k) / dk a ) 2 . This term occurred due to the 
fact that the Peierls substitution was made in the initial Hamiltonians ([5]) and (fl"U|) 
rather than after the diagonalization of Eq. (JTUJ) . The second comment on Eq. (|124j) 
is that (t qq ) vanishes if e(k) is taken in the linear approximation. This reflects the 
absence of the diamagnetic term in the Dirac approximation. The co rrect way is to 
firstly take the derivatives in Eq. (|124p . This leads to the final resullP^] 

¥■ - L w? K(e(k)) - nFhe{m e(k) - (125) 

Eq. (|125p is equivalent to Eq. (|124|) . Note that (t qq ) is always positive and does 
not depend on the arbitrary choice of the sign before t in Eq. ([5]). Now Eq. (| 1 25[) 
is e 2 /h 2 times — 2/(3v3)('~ —0.39) of the kinetic energy per atom instead of —1/2 
for the usual square lattice. 

It is useful to separate explicitly the contribution (t xx (/j, = T = 0)) of the Dirac 
sea from Eq. Ijl25j) : 



where 



(t xx ) = {t xx {h = T = 0)) + (7^0*,T)), (126) 
(t xx (p = T = 0)) _ e 2 a 2 f d 2 k 



V 3ft 2 J BZ (2tt) 



Y2 



(-e(k)) (127) 
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is the contribution of the Dirac sea (the energy of the filled valence band) and 
ft£(/i,T)) eV f d 2 k 



r ; ( , );)2 K(e(k)) + l-n F (-e(k))]e(k) (128) 

is the electron-hole contribution. The numerical calculation of the Dirac sea contri- 
bution (|127p with the full dispersion (TT2|) gives 

(r^ = r = Q)) =a g QRjQ61 

V nr 

The same answer also follows from the linearized Dirac approximation with the 
trigonal density of states if the band width W is given by Eq. (|4"2")) . The electron- 
hole contribution (|128p can be estimated analytically in the linear approximation 
for the dispersion law. In particular, for fi — the temperature dependence of 

the diagon al conductivity sum rule is ~ T 3 , in contrast to the often found T 2 
dependence^llimmimi 

in tight-binding model of conventional quasiparticles. 
On the other hand, for |//| ^> T one obtains 



V ' 97r^ 4 w 2 



[| M | 3 + ^ 2 HT 2 ]. (130) 



We note that the \p\ 3 behavior can be observed in an experimental configuration 
which has very recently been used in Ref. [2j[ by incorporating the specimen into a 
field-effect device. In this case the chemical potential p, is easily changed by varying 
the gate voltage V g . Unless /i is chosen to be large, the change in the sum rule 
from its p = value is small [~ (|/x|/£) 3 ]. Using Eq. (1132p for carrier imbalance 
p which is proportional to V g , we obtain that p oc sgn V g •\/T^l and, therefore, 

(T? x (v g ))/v~-\v g \*/\ 

4.4.3. Hall conductivity sum rule 

The RHS of the Hall-angle sum rule (|123[) for graphen o^^ j s equal to 

1 eBta 2 2 , 

Wff = -S7? pfl - (131) 

Here p is the carrier imbalance (p = n e — nu), where n e and rih are the densities 
of electrons and holes, a is the numerical constant from Eq. Q129p . The carrier 
imbalance for B = T = and in the absence of impurities is 

p = ^¥vJ- (132) 

Since ta 2 /h 2 has the dimensionality of the inverse mass and pa 2 is dimension- 
less, Eq. (I13ip has the correct dimensionality of the Hall frequency. Substitut- 
ing Eq. (|132p . expressing vf via t and employing the Landau scale L(B) used in 
Eq. (|108p . one can rewrite 

. H = - A -^lL 2 (B)^. (133) 
yira hi- 3 
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As mentioned in Sec. I4.4.H in the recent interpretation of Shubnikov de Haas 
measurements, a gate voltage-dependent cyclotron mass was introduced I 14 | 15 | 
through the relationship \fi\ = m c v F . If this is used in Eq. (|133p . we get 



eB 

u H = 



Since a full upper Dirac band corresponds to a value fi — W — \J y/Znt (see 
Eq. (22)), in this case formula (|134|) resembles the formula ujh = u c = eB/mc 
from Ref. 1 1131 for a two-dimensional electron gas with m c replacing the free electron 
mass. In graphene, however, m c varies as the square root of the carrier imbalance 
\p | and the two cases look the same only formally. 

When the spectrum becomes gapped with E — ±^Jh 2 v F \i 2 + A 2 , where A is 
the magnitude of one of the Dirac masses (|9ip . the carrier imbalance is 

(/i 2 - A 2 )<% 2 - A 2 )sgn^. (135) 



7rh 2 Vp 



This implies that the gap A can be extracted from the change in u>h obtained from 
magneto-optical measurement s. Th is kind of measurement which reveals gapped 
behavior has already been 

donjnsi 

on the underdoped high-temperature supercon- 
ductor YBa 2 Cu 3 06+x- 



5. Conclusion 

The number of articles devoted to graphene grows v ery fast and the seminal exper- 
imental papers where the IQHE was reportecPHS have already been cited almost 
300 times in arXiv, while the total list of references of our review has about 100 
papers. This is partly caused by the restriction on the size of the manuscript and 
by the amount of time we decided to invest in this project. We intentionally omit- 
ted such interesting and important topics as bilayer and multilayer graphene, the 
role of impurities, defects and their gauge field description, minimal conductivity at 
the Dirac point, edge states, weak (anti)localization, fractional quantum Hall effect, 
graphene nanoribbons, etc. Each of these topics may, in the future, become a sub- 
ject of a review on its own. Instead we decided to restrict ourselves to the specific 
topics covered in the review and to discuss them in detail. We hope that this cho ice 
is indeed complementary to already existing textbooks and reviews 

During our work on the review, we found that in the literature on graphene, the 
word "chirality" is becoming as fashionable as "Berry's phase", but depending on 
researcher's background, it is used to denote different concepts, so that it is getting 
hard to grasp its meaning. Here we followed the commonly accepted definition in 
field theory-^ and demonstrated that the chirality quantum number corresponds 
to the valley index in graphene. In Sec. 13.41 we also introduced the pseudohelicity 
which characterizes the projection of a quasiparticle pseudospin on the direction of 
its momentum and explained it spccificifically in 2+1 dimension and the relationship 



February 1, _ 2008 11:55 WSPC/INSTRUCTION 

review-graphene'jversion 



FILE 



GRAPHENE: FROM TIGHT-BINDING MODEL TO QED 2+ i 41 

between pseudohelicity and chirality. We hope that this solid state reincarnation of 
chirality will lead, in future, to a deeper understanding of chirality in field theory. 

It is also very instructive to see how the spinors emerge from the two atom 
per unit cell description of graphene. Because these spinors are not related to a 
real spin, the corresponding degree of freedom is called pseudospin. This solid state 
realization of pseudospin may help us understand the origin of a real spin better. 

Also we mentioned a link between the QED2+1 description of graphene and d- 
wave superconductivity. In both cases, the low-energy quasiparticle excitations are 
described by massless QED2+1. Massless excitations in graphene are often compared 
with massless neutrinos. However, as we pointed out in Sec. I3.5.4( this analogy in 
not complete, because their discrete symmetries are different. Massless QED2+1 
has U (4) symmetry and a rather common way in the field theory to break this 
symmetry is a dynamical generation of the Dirac mass. It is unclear at this stage 
whether this rather generic mechanism works in graphene or the observed new ±1 
QHE states can be explained by other mechanisms. From our point of view, this 
is one of the most interesting open questions which demands further experimental 
and theoretical work. Accordingly, in our review we simply tried to facilitate the 
future theoretical work by providing a classification for possible Dirac masses. 
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Appendix A. Solution of the Dirac equation in the Landau gauge 

The Dirac equation in symmetric gauge and in the representation (|52p is solved, 
for example, in Appendix D of the second paper in Ref . [TT] (see also Ref. 11171 for a 
solution of 3D Dirac equation in the Landau and symmetric gauges). For eB, A > 
in the Landau gauge A = (0, Bx) the positive (tp^) an d negative (ip^ ■*) energy 
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solutions of Eq. (|104|) at K + point are 



^\ X ,y;n, P ) = * e -i\Wr»-L= ( VW^»-f I „ > 1, 



4>, y ;n = 0,p) = (_.^) , n= o. (A.1) 



The solutions of Eq. (|105[) at K_ point are 

1 c -i\E \t+i P y I 



4 +) ( a; ,y;ri = ! p) = ^= e - i l^l t+ ™f . ° ), n = 0, (A.2) 



n > 1, 



Here £ = x/l + pi, where I = (hc/leB]) 1 / 2 is the magnetic length and p is the 
momentum along the y axis, 



E n = ±J A 2 + hv 2 F 2n\eB\/c, (A.3) 
«>»(£) = (^ 1/2 2"n!)- 1 / 2 e -« 2 / 2 i/„(0, / dtw n (Z)w m (£) = S nm 



with H n (£) being the Hermite polynomial and we defined = 0. The lowest 

Landau level n = is special: while at n > 1, there are solutions corresponding 
to both electron (E n > 0) and hole (E n < 0) states at both K + and K points, 
the solution with n — describes holes at K + and electrons at K point. The 
corresponding structure of Landau levels is shown in Fig. [3] (c). See also Fig. 2 of 
Ref. 11181 for a graphical representation of the dependence of the LLL position on 
the relative signs of eB and A. One can also say that for eB > the solution 
of Eq. (|104p for LLL at K+ describes holes that live on the B sublattice, while 
the solution of Eq. (|105[) for LLL at K_ corresponds to electrons living on the A 
sublattice (we remind that in the K_ valley, the spinor components are inverted in 
comparison to the spinor in the K + valley). Finally, we notice that the relation in 
Eq. (|34|) is violated in the finite field. Indeed 



T 2 n* {k,eB)T2 = -HQ(k,~eB), H (k, eB) = Tik x + T 2 (k v + eBx) + Ar 3 . (A.4) 

This indicates that the wave function Ta|^>)* is the solution of the equation 
Ho(k,-eB)T 2 \ip)* = -Et 2 \i{j)*, i.e. the generalization (|A.4[) of the relation (f3"I)l 
for a finite magnetic field does not describe the symmetry of the spectrum about 
E = 0, but rather relates the solutions of Eqs. (|104|) and (|105[) obtained for eB > 
with the solutions of the same equations for eB < 0. 
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